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Preface 



The present booklet describes fourteen elementary experiments in Mechanical Engineering, of < which 
more than half can be readily Incorporated into a sophomore course. The book di fferrr signif icantly from 
most laboratory manuals in its philosophy towards' both the student and the instructor. 

The student is introduced to the concepts of experimentation, to the inter-relation of theory and 

« 

experiment, to the use of dimensionless quantities and to the vagaries of instrumentation with a minimum 

/ 

of prior theoretical background. Some of the experiments are open-ended oo that students can learn 
early in their college career that a- slight change in operating conditions can shift a problem from the 
trivial to the challenging. Perhaps most important is that many of the experiments involve transient 
situations; all too rarely are students* introduced to the idea that the real world does not correspond 
to the steady state ^bich is usually ,set up in the laboratory. - • 

The instructor should appreciate the underlying simplicity of the experiments; little is demanded 
in terms, of manufacturing and erection, and the instruments called for should be available in any 
undergraduate -laboratory. Where needed, construction details are given in full. Furthermore sample 
results are presented; all experiments have been pretested with student groups. 

In terms of required background knowledge, little is assumed for experiments 6, 7 and 12 which can 
be offered at any time during the sophomore year and may even be included in a freshman engineering 
course. Some knowledge of differential equations is desirable for experiments 1 through 5 with, a dyna- 
mics course recommended for 2 .and 3 and strength of materials for I and h. Experiments 8 through 11, 
13 and Ik are probably' best offered during the junior year although for exceptional students or for 
special .projects they might be suitable on the sophomore level. • , 

Hopefully this" set of experiments will make an impact in a field where all too long stereotyped 
performance tests have ruled the curriculum and where experimentation has been driven from the lower 
division years.' If it can provide only a modicum of excitement about mechanical engineering problems 

i 

and supplement the analytical diet of which so many sophomore complain, it will have served its purpose. 

- ^Th is booklet has been prepared under the sponsorship of the National Science Foundation 
(Grant No. G^2h6f) and the continued assistance and patience of the Foundati6n is grateAilly acknowledg 
One free copy is being -distributed to 1 each accredited Mechanical Engineering Department. Additional 
ctfpie;s will be. available at w nominal, cost from the project director. 

; " % ' > F. Landi6, Chairman 

. " ■ , . " , Mechanical Engineering Department 

- • . ^ New York University 

February 1971 V . \ Project Director * 



Introduction 



The objective of tKis work i& to bring together a series of experiments in mechanical 
engineering which are understandable to students in their first and second years of the 
undergraduate program. While an effort has been made to maintain this leve^it is to be 
expected that some sections will be too advanced for freshmtn and sophomores- This is most 
apt to pccur £o the theoretical development , and in sections where one might expect this to 
happen there are notations in the text. Each experiment is relatively simplfe to construct 
arid operate, and cfcn be done using instrumentation normally available in undergraduate lab- 
oratories. Therve is no attempt to prepare a structured laboratory program- The order of 
periments is arbitrary except in so far as they are grouped aC^^ding to fields of study * 
^ Each experiment is written in sufficient detail to be useful to the instructoryor 'the 
stu Hk nt - Primarily, the text is designed as an instructor's manual. ' Theory is developed 
with ^ftpnciTJerable detail when that is necessary, and usually references are givqn for two 
sourcel^where additional pertinent information can be- found . Drawings are given of special 
Gfluipmeffit as well as diagrams for the experimental set-ups. 

Eaciy experiment hT^s been executed successfully by 'undergraduates in the Department of 
Mechanical Engineering \jat Tufts University. Sample results, curves, data and calculations 
are given. D.ata can not always be guaranteed to be precise, but it is safe to s^y it is 
the best available from several' sets of results. An object statement -precedes each experi- 
ment. An effort has been made to insure that each description will stand by itself. It i$ 

i 

not necessary to read t\\e experiments that precede a particular experiment in order to under- 
stand it. v 

An objective of any experiment for the .engineering undergraduate is to motivate the 
student ;and call On him to. use his imagination and to be creative. All of the experin\ents 
can be modified so that the same apparatus can be used to do several distinct^experiments 
or a Succession df experiments. It would be improper to use the same experiment repeatedly 
for a succession of students, and that plainly is not the intent of setting these down. 
Some modifications are given in each of the experiments. A useful technique for introducing 
modifications is to have one group of students start an experiment, letting groups that 
follow build on the results of those who precede. In this way, the development of the experi 
ment follows .the pattern of long-term research programs. It also, provides incentive for each 
group to present its results in a manner easily understood by the succeeding grotip. In the 
author's opinion, these experiment's are best done with three students in a group-. 

In addition to the people cited In the text, I sho ill* like tcf Acknowledge with appreci- 
ation the speoial contributions of several: Professor Fred Landis, for his confidence in me 
which he showed by asking me to undertake this task; Peter Boyce, v> Tufts undergraduate, for 
preparing the drawings; Mrs. Irma Wallace and Mrs. Betty Steel for typing drafts and final 
copy; and to the students too numerous to list who performed the experiments so well/ 
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Experiment No/ 1 r " 

Tho Tors ional Pendulu m 

^ * 

The object of this experiment is to - measure the frequency 
A Of a torsional pendulum and ^determine the stiffness of the 

system. * 

A simple, but useful piece of laboratory equipment-, is a length of music wire about six, 
feet long with screw fasteners brazed on each end, Fig. lAa. One end is clamped in the 
ceiling, while the other is fastened to an object to form a torsional pendulum . A torsional 
pendulum oscillates about an axis, in this case the V-ire , as shown in Fig. 1-jlb. Several^ 
interesting experiments can be developed from this- system, but first the wir4 stiffness must 
be determined-.- -, - • — -' - - 



.Wire calibration , v 

To be useful, the stiffness of the wire must be known. s For kmall oscillations, with / 
negligible air resistance, the motion of a torsional pendulum can be described by the equation 

ili fko = o' , ' V 1-1 

dt 



where I = mass moment of the object about the axis 

of rotation 



9 = the angle of oscillation 
t - t ime > * 

■»> 

-k = the stiffness "of the wire, radians per unit: 
of torque - ^ 

From the solution of this equation which is 9 = A sin (u^t + a) it is seen that the dis- 
placement is sinusoidal with a frequency, called the natural frequency, of u> n -jfy ..radians/ 
second. (See'ref. 1). In this 'expression , A is the amplitude and ex a phase angle, found 
using the initial conditions. Using this information it is apparent that if we measure *he 
frequency of Jie oscillation for an. object of known inertia, we can calculate -<he stiffness k 
A calibration disc, Fig. 1-lc , was made from a solid cylinder of metal (aluminum). Its 
mass moment of inertia about its axis of rotation^ is 

• • T - 1-2 

• t . \ 2 • 

where m is the mass Jf the disc and R rts radius^ Measuring thp period of oscillation of 
the pendulum, the stiffness of the wire, k, can be determined from the expression for t^. 
There is a discussion of units included with sample results in the next section. ^ 

An assumption has been made that the mass-moment of inertia of the wire is. negligible 
when compared to the inertia of the calibration disc. Each section of the wire turns through 
an angle between 0 and 9 depending on its distance from the fastening -point Hv the ceil- 
in*! Normally an "equivalent inertia" is defined as the inertia which would- be added to the 
disc to produce the same effect as the inertia of the wire. As the wire turns through an 
angle less than 9, this is bound to be smaller than the inertia of the wire, l wlre - ,}!hl - s .. 
equivalent inertia could be found experimentally. If, instead of I in eq. (1-1), wo add 



--*| 1/2 
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1/4-20 th'd: 
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the equivalent 'inertia to that of the' calibration disc, I, then 



= \l k 1 



-co ik i-3 



Now <can be measured, I "»an be computed f ro m eq. 1-2, leaving two unknowns, I q and k 

Using a second disc of known inertia 1^ and repeating the measurement for^^, produces a 
second equation in k and 1 , which can be solved simultaneously with eq, (1-3). It is 
possible that" the effect of 1^ on the frequency is too small to measure accurately, as is 
quite likely the case here. It is a simple- matter to determine I analytically for a wire 



or shaft of circular cross section (reference 2).- The result is 

I = wire / 1-4 



3. 



An assumption has been made in this analysis that the angular deflection of the wirfc is pro- 
portional,!/^ its distance from the point of . attachment . 

If as possible to compute k from the properties of the music wire using the shear 
modulus for the wire G = s/T, in lbf/m , . . « 

where <j s = shear stress and T = shear strain. 

Let ( = length of wire 

0 = angle of twist per unit length of 

A = c^oss-sectional area of the wire 

r = radius *of wire \ 

then the moment required to twist a unit lengtfi of wire through an angle 0 is 



/ r sdA 

A 



/ 



r- = J v GYdA 
' A 

y = r0 

S 2 
Hence M = I " flr GdA 



i 



r i 

M = j»Gj" 

Where J is the polar moment oi inertia of the wire cross section 

Now k = 7T t torque per angle of twist , where 9 is the angle 

9 



"through which tKe end of the wire 'is twisted. In terms of 9, the deflection of each 
*#Rbment is 0 = Q/l 

0GJ * GJ i- 5 



Hence 



£9 f. 

Sample results • , ■ . 

An experiment 'was conducted to determine the wire stiffness. 



- 3 



Dlmanaiorus of calibration* disc: , 

* ., mass; -m • 1671 r gm or 3. $01" lbm 

^ radius:- R ^ -9 cm ± 0.1 cm or 0.295 : ftJ 
, ' ^ I=l/2mR 2 > = } l/2-(1671)(0)?" * 6.77 x 10* 4 gm-cm 2 • 
I b (l/2)O.68)(0.295) 2 = 0.16Q ibm-ft 2 . 



ok, in English ui\its 

' v For -five measurements of the period made with a stop watch, the average was 

T 22 18. M sec. which corresponds to frequency of 

a) = — s 6.3 l » radians/Second 
n T f i 

m *■ ' ■* » 

k ^ a) 2 I = (0".3M) 2 (6.7?-10 M ) • 

n - , ; ; - 

" - M 

= 0.783 x 10 dyne-om/rad 

k 



or 



and 



( ; 3 i4) 2 i2J^i(i2) B 0.00609 lbf-in/rad; 



32. % 

This calo^ation required* introduction of the two 'conversion factory / 

. 2 

J - -g =? 32v2 lbip-ft/ibf-sec . >' . - ^ 

.12 in/ft to produce the units ordinarily \i^ed for k, 
A value of k c^n "be" found using- eq, (1-5). For the wire: 

m length : 



diameter 



# -± 1M8 cm±0.1 cm -or 
* ' 58.27 in ±0".0H.«In 

d « .6* 10^.;lo" 2 dm ^ , 
2.M -x 10-? in 



J = ird 



32 



r Il(6/l0)''l(^ 8 

3*2 



om • . op 



\ 



(| - b 
lt(2.«0 10 ..4 
52 1,1 



Shear modulus for music wi 



\ 



From eq. (Ir5) 
In cgs: 



Q - 



k = 



' r o 
12.5 x 10 lbf/in of 




t . 2 



lbf/cm )(»K>t5 x 10,.dynes/lbf)<0.155 in /cm ) 



. cm 



Gj. ' ( 8.62 x 10")ir(Q*.Q) 10 
'St - . (iM-8)(32) 

0.792 x 10 dyne -cm/ rad 



In English units : 



k «= 



(12.5 x 10 6 >TT(2'.M) 4 10 Q 



(59)(32) 

= 0.00690 'Ubf-in/rad . 

The experimental value is within 1% of these predictions. 

To'.deterroine the effect of the mass-moment^of inertia of tfle wire first compute 

♦ ' (m . Mr : ) 

wire wire 



2 



. For thib we require the 



tt(^:u) 2 10^ M 59(5001 



Mass df the wire: m , * ■■ (— O(density) - — (l72B) 

7 " r^f .' * ■ j - 

3 7.72 xl(T^ lbm I 

/ 



.,The density cho?en,for the stae f Wwire was 500 lbrf/ft . . . 

• ^ t (7,72)(1.2)^0' >t . Vj ^ , , 



Try 



jr. . . wire 

' . , ' ;■ • ' 

.' ' 

'Then determine the effect" on, the inertia by ; introducing. 



• i 



'V Iwif ^ = '1.85 ^lO"" 



t V I ' = V 0.160 + .0002 " 07,1602 



....... \ 



!- p The co^ected inertia inc luring this is v 

.- which -is'a change in the value of k of about 0.1%: ;This la small enough to neglect. 
, • » References 



1. Church, A. 1, Mechanical Vibrati ons, John Wiley and Sons, Inc., New York 
. <\ (1963) p23. " ' 

2. Crossley, F.R.E- - D ynamics in Machines , The Ronald Press, New York (195M), p!19 
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* .* * Experiment No. 2 



Damping, a Torsional ltenflulum " * 

• , An Electric galvanometer consists of a cylindrical oore suspended on a filament. When ^ 
the core moves, it carries an indicator such as a needle or mirror.- It is necessary Jo damp 
the system so that it does not oscillate when %ie reading is changed, but comes to rest at 
the proper position . . One suggestion is to immerse pdrt of the core in a fluid so that the 
shear force between the fluid and the wails' of the cylinder will provide damping . 

1 ■ « * 

The object of this experiment is to obtain a relation between v 

the geometry "and the damping, as ,-we.U as' between fluid properties * ^ 

• * and dampings In addition, a second experiment will be- proposed . - 

to Qxami.ue_the effect of a stationary concentric cylindrioal_wall 

on the doping. . „ * 

Experiment No. 2 A 

t 

Effect of Fluid Properties 

' 'A model' of the system Is shown in Fig. - It consists of a hollow cylinder, 3.867 in. 

diameter/ M ll/e'/'in. long" A It is capped with a threaded hole to receive the end of the music 
wire. The .cylinder weight is 2.321 Ibm. A drawing with dimensions is given in Fig. ,2-2. A 
paper is cemented on the top surface, having markings on the outer edge in degrees. This will 
be used to measure the'tengle of oscillation, by matching the lines with a stationary indicator 
at- each .end of the oscillation. 

* -Using the music wire, the cylinder is allowed to oscillate in air and its period of . 
oscillation determined. From the. method described in Experiment 1 (which assumes air damping 
to be negligible) the inertia of the cylinder is found.* ■ ^ ' . • ' 

-The .pendulum is then allowed to oscillate in a liquid, as illustrated in Fig.- 2-1, and 
.the angular displacement-measured as a function of time. For the fixed geometry of the cyl- 
inder, the variables of interest are: the kind) of fluid used; the fepth of immersion; and the 
distance .between the' cylinder wall and the container wall. In the practical case it was the 
last ,ef feet thatrwap <$$ npst interest. 



* ™ > • - . 

Deve%opmerit of t heory-. - • \ ' 

, ■ ^ ... ■ 

The. objective was^to ,de,$jarmine the damping ratio for the damped single-degree system, 
j.fhis can be' expressed of the .logarithmic decrement, .5, or the ratio of "damped fre- 

^ency/ W|V Jto the-^natural frequency, u> n . Damping ratio is the ratio between the coeffic- 
ient 6f v.fecous damping, c, and thl critical damping coefficient c^,.. The development of * 
tifcse expressions can be found in any elementary book on vibrations, e.g., ret. i<or /. 

A To summarize this analysis briefly, begin with a statement of 'Newton's second law for 
\the system'- Equaling the time rate of change of angular momentum (16) to the sum of the 
' resist ingfeiip due to the spring (~k6) and the damping (-c6) results in .the governing 
' differential, equation 

«' • ta . _A i - n 2-1 



3.067 
DIA. 



FIG. 2-2 HOLLOW CYLINDER 



FIG. 2-1 EXPERIMENTAL MODEL OF 
GALVANOMETER DAMPING DEVICE. 
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FIG. 2-3 ENVELOPE OR . 
AMPLITUDE DECREASE 
IN AN UN PER DAM PEP FLUID 
IN TWO DIFFERENT CONTAINERS 

1500 ml CpNTAINER C/C c * 0.0338 
3000ml CONTAINER C/C c = O 0339 



4 — »-TIME 



is 



wher. the dotted 0»s are 4 the second and first derivative* with reepaot to tif. Introducing ^ . ^ 
o "trial solution of the form *' , x 



leads to th« "solution 



where 
and 



0 ■ 


A P 1 

V 


0 " 


V * 

A.« P 1 






h'- 






G 







s. 



.Pot 



2-2 



2-3 




2-fr\ . , : 

• v. . i 



c .2 k 

•IF/"' i. v ^ 

The 'form of the solution depends on \:hV hatbre of . the term in the radical, ..Wienjt, ' / 

\- t • ■•■ - ' -I ' 7 • ■ < A , 

aero-, -• . " ■ - ... .Y-fra -.: . * . v ..... . . : _..*'..:../. . ... !. ..: 

' " " * ' T - .-/- ./ 

The system is said to be critically damped and" the critical damping factor if J 

- . ; | '' 

, c - 2 kl . h 



/ 



;2-5" 



2-6 



This can be introduced Into the expressions for p , by multiplying and divijfrind by 



The result is 



P 



1.2 



| 2Vkl f M kl • Jl 



2-8' 



2_ ± 

c 

c 



The ratio c/c is soifietimes expressed as C 

When c/c >1, the system is said to be overda.nped , and .does not oscillatL Instead 
the amplitude decays exponentially froln the initial -state. ' Light damping or unWamping 
occurs when c/c^l. Then the root is imaginary and the solution has the form c|f a damped 
sine wave . 



(-)«t 



to 



0 = 0 



o to 



sin^l - (f-) 2 w n t + *] 



Here 0 is the init ial' amplitude which decays in time according to the exponential exprqs- 
o & 



- )to t 



sion e 



c c - n * 



The frequency is called the damped frequency, w d ,_ and is 



03 



-V 



1 - eo 



2-9 1 



Aefined to be 
2-10 

, the phase -. 



This is the actual frequency of the damped system, and can be .measured . - Finallfc 
angle $ can be determined from initial conditions. » 

Solving the expression for the damped frequency provides a means of determining th^ 



ERJC 



6 



If w n is frhfc frequency of the^pantfulum- in air (undamped) and the mqasiu*fcd frequency 

in the damping' .fluid , then *.c/c can be computed. 

Whale th* ^preceding procedure 'can 'be "smploybd where "damping .is significant it is 
di^fioult^to determine exper£m£ntfi(lly when damping. is smAll because is approxi- ^ 

matoay' equal to* , Aft alternate pfoc&duVe is to measure thtf amplitude -of^tha oscillation,, 
0, on two succt\3yiv6. Oscillations, 0 and ' 0 ? . - L*tf: t>« swings l^e at tiroes t^ J*™ 

t v + 1 where ^ T ifc the period of the damped -penduium. 'Then since the sine pjprt of the 

1 _ / » _« , 1 1 2* 

.expression is the same i^/edfeh case, the ratio 6F the two amplitudes' Is 



<57 



c 

— 0) t, 

c n 1 

c 



CO T 



— 0), (t- tx) 
c n 1 
c 



2-12 



The logarithmic decrement 6 is defined as the logai^hm of the ratib of successive ampli- 
tudes, or 



6 = JLn 



- 8 1 



2-13 



Tor a more accurate experimental determination of 6, it is better to compare two peak, 

7 • fa 

amplitudes that are several oscillates apart, say n . Then it is easy 1o show that the 
logarithmic decrement between the peak and the (i + n)th peajt becomes 



i+n 



2-1H 



But the period is relate/ tp the damped frequency by the expression 



2ji 
w d 



2-15 



With ttyis substituted in th^' expression for th<? logarithmic decrement, the damping ratio can 
be expressed 'as 



c = 
c 

c 



± 



2-16 



* Sample results 

For the cyii/der oscillating in air, the damping was too small to determine the ampli- 
tude change, so pie value of 6 was assumed to be small and a> set as a> n . For the system 

tested a> = S^,M rad/mini The inertia of the cylinder was found to be U , M6 Ibm-in , using 
n / • i . 

^the music wire^ calibrated in Experiment 1. 

In one/experiment water and clycerin were compared. »lt was found that for^water, the 

damping was^small and it was more accurate to measure amplitude, while changes in frequency 

of the glycerin-damped system were more accurately obtained. Some results are given in the / 



following table, for tests Un containers of several sixes. The effect orthe container wall 
is evident An tti& 'results. 



( 




Fluid • m eontftir^ V - * n . 2> ' * . '- c 



Wat^r ' ^ 5 pfciir ; ' *V v .2Jk3^' ^; \\ ■' . '.QSMQ^ 

- V; ' 1500' ml beaker' 1 % \ , * $H .•• x - .03^8 , 

s- " . 3000 miybeak^ -j 1*7/ ;/ - < . -.9339 

Glycerin 5' qt . pail. V^'"' 31-9 rad/m,in*. .037 

1500 ml beaker 20.2 rad/min • .938, . 



■ -The glyeevin is a mo^ off active damping fluid than water and the effect of the^w^l ia 
\ndre pronounced. T>ie^maller gap (.1500 mi beaker) produced more effective damping than the 
Jarger (5 qt . pail). Curves of amplitude and the envelope of the "amplitude are shown In 
Fig. 2-3 for the 1!>00 ml beaker. - v * ' 

Experiment 2B ' 

■ 

Effect of a Concentri c- Ha ll » „ 
. 1 — $ 

In the proposed galvanometer system, a concentric wall will be employed to improve the 
damping. The object of this experiment is to determine the relation between damping and wall 
location. The cylinder is allowed to oscillate in glycerin in a large container. A concen- 
tric wall is introduced as shown in Fig. 2-Lf In this case the wall can l?e made from a. sheet 
of thin aluminum so it can be adjusted to 'several diameters. -As the cylinder is 3,867 in. In 
diameter, the gap Will be half the difference between wall diameter imd 3.067. Tests- then can 
be made by determining amplitude a& a*function of* time for the torsional oscillation, for 
several V/all diameters . \ ■ ♦ 

* ! 

So me resul ts " > . 

Damping' in, these tests was more severe than in those for Experiment 2Ao It was over- 
damped so tMt the pendulum did not oscillates . Amplitude was ready by several observers at 
.intervals and recorded. A more accurate measure wouid be to use stroboscopic light and 
record the position by multiple exposures on camera film. The angular displacements could 
then be 'taken from photographs. - l i. *, ,. v , * . m - 

Curves of' displacement versus • time "( IMg . v 2 -M ) show the interes ting result of an optimum 
wall diameter of 5 inches. Larger diameters lead to a partial oscillation whilo smaller 4 
diameters are so heavily damped \j>at the cylinder requires too much time to c^itie to rest. 
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FIG. 2-4 EFFECT OF A CONCENTRIC WALL ON THE DAMPING OF A TORSIONAL 
PENDULUM, THE FLUID IS GLYCERINE. WALL DIAMETERS ARE SHOWN 
ON THE CURVES. CYLINDER DIAMETER « 3.867 in. 



. - s ' , Experiment ^o. 3 

t 

jU j^ffect of riu id_J^reeB_onJta^^ Systems 

the object of -this experiment' is to investigate the effect a 
sloshing fluid hws 4>n a dynamic system. 



Introduction 



In the previous experiments the pendulum was a Solid . mass , and the position of each . 
element in the mass was fixed with" respect to each other. If, during the oscillation, the 
pendulum interacts. with a fluid-, the fluid will be sej in motion through viscous interaction 
and pressure force's- .The fluid will then move '.at a speed which is different from the speed 
of the pend^um. ^In solving equations of motion, the mass of . the fluid , so excited, tfll 

-increase the .mass \>r- inertia term as well as contribute to the damping term. -Rut 'as the 

velocities of. the fluid and solid surface are different, the effective increase i<f inertia 
will be different from the actual mass or inertia of thft fluid. It cah be described- as an _ 
equivalent mass or inertia. This is a difficult problem to analyze but experimental results 
will be useful in understanding the behavior. In addition, these ' experiments provide an 
opportunity to examine dimensionless numbers. * • \ 

Dimensionless numbers allow generalization from, an experimental result. For example, 
if one says a measured length is accurate to within one foot, it m*y not carry much meaning, 
as it would be good or poor depending on the length measured. . If, instead, the error was 
one foot in a length Of U000 feet, then one can expri!^ this error as a dimensionless num- 
ber, 1/M000, which gives an indication of the .accuracy . A more complicated expression is the 
drag "coefficient of a body moving through a fluid medium- This is the ratio of the resistance 
force created by the fluid to the force equivalent of the velocity pressure ^~ - Using 
C D to represent the drag coefficient and A, a projected frontal area of the body ,~this is 
traditionally written 4 



r n rag • 3-1 

D ~ 2 

(^) A 



In dimensions of force (O , mass (M), length (L) and time -<T) we see that 
x Drag , =" force , F 



and 



3 

p = density of the fluid, M/L . 

v = velocity, L/T 
■2 

A = area , L. 



C n ■ F . 2 _ fV 3 ._2 

D = _ ML 



(\> jk) 



There is a redundancy here as we know. that- force , mass, length and time are related- by Newton's 




- FIG 3-1 A TORSIONAL PENDULUM USING A BUCKET , 
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Consequently 



n,-2 



'D " T 2 ML' 



c ^ B !% (L_) = dimensionless. 

v 

Reynolds number 



pvD . - 3-3 

w ith ^ D • = a length (L) 

V» = viscosity ML/T 

is another useful dimensionless expression in fluid ^chanics. Experimental results fof the . 
fluid dnag on bodies /such as 'a sphere , show that a" curve, of versus Re is independent ' 
of the size of the sphere or the properties of the fluid. Thus results from a single experi- 
merit can be applied to^pther similar -systems. • ----- - 

Experiment 3A ^ . f A ^ 

Fluid Mechanics in a Torsional Bendulum* 

A simple but effective experiment can be made using a common bucket as a torsional pen- 
dulum v The handles of the bucket are brazed to prevent swinging, and a clamp of some sort 
fastened to the handle. Fitting the .clamp with a 1/4-20 thread allows it to be attached to 
the music' wire of -Experiment 1 to form a torsional pendulum.. tW. is illustrated in Fig. 3-1. 
'Having calibrated the stiffness of the wire, the mass-moment of inertia of the pail can be 
determined using the procedure of Experiment 1. t 

Next, the bucket can be filled to different levels with sand or other granular material. 
The inertia can be determined using the music wire, and might be graphed against the maa.s of 
the sand in the bucket- It is possible to separate the contribution to the total ' inertia due 
to the sand. It would be best if a granular substance having the density of water were used 
to compare these results to the next part of the experiment using water. 

A more interesting experiment results when the bucket is filled with liquid. This prob- 
lem is a difficult one,. but it offers an opportunity for the student to select dimensionless 
variables and determine their appropriateness. One can define a linear damping • factor , c, 
from eq. (2-1). Using inertia determined from the sand experiments, and, the wire stiffness 
evaluate c as in Experiment 2. It is found from the experiment that this is very nearly 
constant for any given run. The damping factor can be related to a characteristic Reynolds 
number and compared graphically to determine if this is a suitable scaling parameter for the 
damping factor. Variations in Reynolds number can be obtained by using fluids of different 
viscosity or density and by filling the bucket 'to different depths. 

A Reynolds number in this case can be defined as the ratio of the torque at the wall due 
to the fluid momentum divided by the wall torque due to viscous forces at the wall. 
Then 

(mass of fluid) (velocity aP-wall) (wall diameter) 
Re ~ (viscous stress) (wall area) (wall diameter) 

*This experiment was based on one by Professor William C. Reynolds, of Stanford University , 
which appears in ref . 1. 




15 - 

9 



3 



Introducing the properties from this -problem the Reynolds number is 

2- m 



D s= bucket dis^neter-J. 

U a viscosity of : liquid, o ' 1 . ^ 

p * w density of liquid,- 

6 = angular velocity of bucket 




As we require a single value /for 0, we can recognize that* - * 

■ / x ' * ' ' i . ' • '■" 

9 - Acosin cot , • l# " - 3-U 

.where A is an ^amplitude and U) the frequency of the oscillation . 

. " t - * - 

R e - r or £ , ... 3_5 

Whether or not the initial displacement of amplitudc^should enter is a question to^be examined 

experimentally. 1 ^ 11 

This experiment provides a good example of a very difficult problem in which dimensional 

0 ■ 
analysis is a profitable tool to employ/ In fact, by experimental measurement, the rather long 

list of important variables can be narrowed. For example, the ratio of spring constant to' 
moment of inertia is important, but not these individual quantities. Fluids of differing 
viscosity and density can be used. In addition, the spring constant can be changed by 'shorten- 
ing^or lengthening the wix-e , and the moment qf inertia changed by adding mass external to the 
bucket, so that the dimensiohless parameters can* be varied in a number of wayS , and the scal- 
ing laws given a thorough test. 

Technique is not too important. With periods of the order of 5-10 seconds it is an easy 

. . - *> 

matter to measure the amplitude at the zero velocity points by marking on polar graph paper, 

and noting the time on a stop watch. One must be careful in locating the clamp to- have it 

ahpve the center of mass* of the bucket, lest lateral swinging motions be excited when the 

f 

torsional pendulum is allowed to oscillate. With water in a one gallon bucket, the amplitude 

is down about 50 percent after about six swings, so the experiments go rapidly. A three-man 

group can get more data t^ian it can process in a matter of an hour or two. 

The flow patterns during the oscillation can be observed by intrQducing dye near the 

> 

wall and along the bottom of the bucket . Through these patters , the relationship between the 
motion of the bucket and the motion of tTie liquid can be examined. 

Experiment 3B 

A _Torsional_ P_ e j^ulum w ith Two Degrees of freedom 

" * ^ ■ i 

It is simple to modify this experiment to one having two degrees of freedom by connect- 
ing \a sjecond bucket similar to the first with a second music wire as shown in Fig. 3-2. This 
is said to be a system of two degrees of freedom because it requires two independent coordin- 
ates to describe the position of elements of the system. In this case, it is the angular 
position of each bucket with respect to* the ground. This system ^as two natural frequencies, 
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whicl) can be observed If 'the -oscillation is started' in just the right way.. Normally t;he 
motion Will J)0 tlio result of combinations of the two "modes" of oscillation Associated with 
the two natural frequencies. Analyses ami examples of freq vibration (no damping) for systems 
of two degrees Of freedom are given in most elementary te^ts in vibration such -as refs. 2 and 3 

Results a&e better- if the two wires have the same or nearly the same stiffness. It will 
.be necessary to determine the stiff notes of both wires by the method of experiment 1 as ^11 -as 
tfje mbment -of inertia of Qach n\ass . .'TJxe- attractive aspect of this system is that; the inertia 
of either mass can t>e changed by simply changing the amouHfc of sand added to, ea,ch bucket. The 
motion will follow the behavior predicted by the vibration analysis quite well ai^l can- be 
^related to the inertia of the two masses and the stiffness of the two wires; that 'Is to .-/the 
natural frequency of each individual ..pendulum , ■ , . ■ 




I 

> / .. 

To. produce a more complex system, damping can be introduced by adding Water to the 

buckets instead .'of sand. This system is too complicated to attempt to analyze but provides 

an opportunity to study - the motion thi'ough experiment . \ 

Experiment 3C , " ^ , . 

fluid Mechanics ill a Simple Pendulum ( 

The object i vc is^to determine the effects on damping of a 
simple, pendulum resulting when the mass of the pendulum 
4 oscillates through a Liquid- 

A design calls' . for damping the mot-ion of -a simple pendulum by allowing the weighted end 
to_osci flate in a liquid. An experimental model was constructed as shown in Fig. 3-3. It 
cons isted of a pendulym ^feci l.lat ing on a hardened steel knife edge taken from an old weighing 
balance. The pcyiduluiir bob was rectangular in cross section. It was constructed so- the pen- . 
dulum bob could be in the position shown or at right angles to it. A large tub was filled 
with water and pl^cced under the pendulum so that the bob was submerged for small angles of 
oscillation. Tes/ts were then, made to measure the f requency and logarithmic decrement In an 
effor>/ to evaluate the ef'fect on. the motion produced by the water. 

1 l*or a simple pendulum with a mass, m, assumed to be concentrated at a length 1 #, 
from the fulcrum, the restoring force due to gravity', g, is* -mg£ sin 0 fi These forces are 
shown in the il lustration , fig. 3- l i. Expressing Newton's law for the pendulum, 

\ni - — — - -mgJl sinO 3-6 
t " 

For small deflections, 0 ^ sinG eq. (3-b) becomes 



. die 
dt 2 -' 1 

for which the natural f roqueney , u) , i*s 
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KNIFE EDGE 
BEARING 



6 375 



STEEL F^ING 



1-68 FEET 




> .ml. • ' 




WEIGHTS 



RING 

ROD 
PLATE 



25 lb. 
53 lb. 
5 75 lb. 



VQLUME OF PLATE 
018 ft 3 
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FIG 3-3 DIAGRAM OF APPA RATUS 
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FIG. 3-4 FORCES IN A SIMPL E 
PENDULUM 



•?5 



20 - 



and iho period 
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The development of this solution directly parallels that for oq.(l-l) (see also rof. M ) , 
This Expression for the period can be compared experimentally to the- period „of the pendulum 
measured in air. Tor. the pendulum of Tig. 3-3, the mass is not concentrated at a point which 
will result in a slightly different period from that given by eq. (3-9). 

When the pendulum oscillates in the water, damping will result due to the resistance of 
the water. In addition, the pendulum will induce some .of the water to move with -.the pendulum, 
which adds to the mass in the equation of motion- The liquid will also affect the restoring 
force. Consequently eq. (2-1) will be written in terms of equivalent mess, damping and 
' restoring force, M c K respectively. Assuming that the equivalent mass is concentrated 
a distance <t from the fulcrum, the moment ot inertia about the fulcrum is M q £ and eq. (2-1) 
becomes, for small 0, 

H i ? 0 + CO + K0 = 0' 3-- 10 
e e e 

Now wo caY express the equivalent mass as the mass of the pendulum M plus' the equivalent * 
mass effect of water, e, 

M. = (M + e) • ' 3-11 

1 G 

. Damping is caused by the resistance of the water to the pendulum motion. This can be 
expressed in terms o.i the drag eoeff .ici en t , which is known for flat plates, see ref.^5 or G. 
Referring back to e{}, (3-1) for. the' definition oi drag coefficient, C n , and recognizing 
that each t'erm In eq.. (3-10) must be a moment about the fulcrum, then 

C pv 2 

C 6 = (Drag) • l = A ■ l 

A. is tlie area of the pendulum normal to the direction of motion. But the velocity of the 
'pendulum bob is 



Hence 



"frotn which We «c on elude that 



£6 



c 6 = c r ; p6 2 A£ 3 
e D s ~ ■ x 



C e - 2 



C„pO£ 3 A 

P_ • » 3-12 



finally the restoring torque involves not only the weight ot*the pendulum but the bouncy _ 
effect due to the liquid. Hence the restoring moment for small 0, 

K e 0 = l" 1 ^ - (Al )pg] 0 
and ' K = [mg - Atpg } l 3 ~ i3 



If is important to verity that, each term in 3J. (3-10) has the same dimensions, 

2-2 ' ■ 

i . e . v ML T 



- n - 



\9 



To relate the quantities expressed above to measurable quantities , can use several 

.expressions from Experiment 2. Defining, for convenience* . 



equation 2^.10 for the damped frequency is 

i 



UK «* Mx - (C) W 



Combining eqs. (2-13), (2-l'0 ,and (2-15), 



* 1 o _ 6 i , 2ttC 
6 a — Jen -it ' o 

n 6 itn . ., VX-iO' 



with eq. (3-10), eq. (3-11) is 



6 = 



U) 2 



Referring now to the expression for c , eq. (2-6) , thd damping factor can be written 

"J 



Noting that U) = a/k 77 
n J e c 



and 



C * c C - 2 V*< I C 
q c fee 



K 



C = 2C — 
e (0 

n 



C . K 
ve e 



5cw 
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The exp eriment- 

The specific objectives of the experiment are to evaluate the equivalent mass of the 
.liquid and the equivalent damping from measurements on the pendulum. Simple physical measure- 
merits were made to determine the pendulum mass, area and volume. The pendulum ^as allowed to 
swing in the water. A Polaroid camera was used and the shutter opened for a swink of the 
pendulum. Aftter several measured swings (n) the camera shutter was opened again. \ The 
exposure in" each case was long enough to ^record the extreme position of the pendulum X This 
showed quite clearly as the arm came, to rest at that point. During this interval, the Ramped 
frequency was measured counting swings with a stop watch. The experiment can be run with\±he 
bob as shown in Fig, 3-3 or a( right angles to it. 

The procedure for determining the equivalent mass and damping were: 

1. With liquid properties determined, the equivalent stiffness, 

K , can be found from eq. (3-13). 
e 

7. Measuring 0 at the two exposures 5 Was calculated from 
eq. (3-lb). 



'30- 

22 - 



Ktiowing 6 and having measured u^, u> "was computed 
from eq. (3- .16) . 

Now C was calculated from eq. *(3-lM) or " (*-l5) as ■ 



f . n 



2 
) < 



5. Having found C, \ and U) n > tha equivalent damping was 
found with eq. (3-17). 

6. With K and w oq." (3-18) was jsolved for I . 

7. From the relation t 

I = (M f e) A 2 
a value for the equivalent mass of the water was found. 



Sample results 



3 

Step 1. For the peudulupi 4hown in Fig. 3-3 oscillating in water (p=62.M lbm/ft ), the 



equivalent stiffness is 



[>,75 - ( .018)(62.M)3 (32^2)(1. 
2 



68) 



250 



lbm-ft 



sec 



(Note- if K is divided by the conversion factor g » 32.2 lbm ft y . . the result is 
e 0 lbf-sec 1 

7.77 lbf-ft, which is properly'' torque . ) 

Step 2. Experiments were made with the pendulum mass as shown in ,Fig. 3-3 and with 
the broad edge of the pendulum normal to the' path of the arc. F*or the latter case, the 
damped frequency**** 25 cycles per minute and the logarithmic decrement measured from the 
1 photographs was 1.20. 

Stejp '3. From equation (3-16), 



u 2 



4> 2 * * 



1 . 101 



"^1. 101-25 * 26.2 cycles/in 
or 2.75 rad/sec. / 



S k tep The damping ratio, equation (3-1$ is* 



V 1.101 



0.303 



Step 5. From (3-17) , the equivalent damping was 



/ • 2 
(2)(0.303)( 250) r e ,* lbm-ft 



or 



5b. 1 
32.2 



1.71 lbf-ft-sec 



r i 
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Step 6, The eqtflUalent inertia, I is from /qV (3-18) 



e 

I * ■ ■ ■ * * ■ ' ■ ? ■ * 33.1 lbm-ft 4 

* K > (*.7sr ' . . . 

whigji i$ equivalent* io a mass of . • 



/ 



(1.60)' 

Stop 7. Mass equivalent of the water is 



• # . a 11,71 - B>75.» 5.9fik lbm 



Based' on the 'results of Stelson .and Mairs (tVf. 7) the equivalent mass, M , should 
have been IS.OJLbm. __Th£s is l'»%_ JLarger than -that S*termijwd .ab£>.ve/„ 

An equivalent drag coefficient can be estimated for the plate using eq.. (3-12) and the 
computed damping coefficient 1.71 lbf-ft-seo.* A velocity is r^uired and th^ average was * 



^assumed to be .w^/yT . * '.^ / 

Then A \ ' / 



(62.ft)(JL.6B)^ 1 ' UM 



= 0.351 

This .is about 1/3 of what' one might anticipate for th><£ drag coefficient from experimental 
results (ref. 5 or 6). This discrepancy could be examined by further experiment . It is 
quite likely that a more appropriate way of selecting, "average velocity" could be made. 

Two additional observations can be made regarding the experiments. Results r for the 
"pendulum turned with its narrow' edge normal to the motion were comparable to the above in 
their agreement with ref. 7, and in the magnitude 'of the discrepancy in determining Cp. 
Several tests were made for each case, varying 8 at the start of the test from 5° to 22°. 
The damped frequency was unaffected for this variation in the starting condition. 
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Experiment No. '* 
Internal Damping in Structural Material 



4/ 

? Prepared'by Robert" Greif 
Associate Professor v of Mechanical Engineering. 
Tufts University 

i 

The object of this experiment io to def ir^analytioally aijd 
measure experimentally the internal damping properties of 
structural materials. * ' 



\ 



Introduct ion ' • • " . 

t This experiment brings together, in a laboratory setting, some concepts in materials 
-and dynamics that -t*e student -learns -duriug-his -first ^ 
Some of the more pertinent aspects are ^ 

i) modeling a beam as an equivalent spring 
14) mgdeling material damping behavior in terms of equivalent 
viscous damping , 

Hi) clear determination of the relative damping properties of 
various common 'engineering materials ar\d the possible 
variation of thi3 damping effect with frequency. 

Nomenclature : • 



m 



- mass of cantilever beam 
M - mass pf tip load 
y - deflection of tip mass 
I - length of beam > 
I - moment of inertia of cross-section 
E - plastic modulus of beam 
t - time 

3 

K equivalent spring constant; <K = 3EI/£ 

c - damping constant , 

- critical damping associated with viscous damping, c q = 2Mw. 



c 



^ undamped natural frequency; u) n = J K/M 



^jlf^ - damped natural frequency; tU) d ~ w n * 

£ - dimensionless damping ratio; C = °/ c c 

n - material damping factor (dimensionless) 

6 - logarithmic decrement' 5 . ■ 

Discussion . 
■ — ' / 

The basic experiment involves the free' transverse vibration of^a tip loaded cantilever 
"beam and the measurement of the resultant amplitude decay. The first step in the analysis is 
to model the cantilever as an equivalent single degree of freedom mass spring system Fig. '1-1. 



■ I -> ..: ... : . , 

M»m 
K-3EI/L 8 

FIG. 4-1 EQUIVALENT SINGLE DEGREE OF FREEDOM SYSTEM 



W *~ T f*rt»t [cos »Jt + ~~ sin o>j 




i 



FIG. 4 -2 RESPONSE OF DAMPED SYSTEM WITH INITIAL DISPLACEMT 
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FIG. 4-3 EXPERIMENTAL SETUP FOR DAMPING EXPERIMENT 
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In the following, it will be assumed that M»m so that the beam is essentially aoting as a 
spring for the tip mass M. 1 The stiffness', or spring constant, K may be derived from 
^•strength of materials" considerations for a cantilever subjected to a tip load. This- load, 
^vided by the deflection at the tip (x = I) due to this load, gives the following value 

of>V (ref. 1) 41 

V K s Clbf/in] . 

\ " ' 

The quantity c is not the usual constant- associated with the viscous damping of a dashp ot 

but is related to the material damping factor n that we are seeking to find in this expert- 

ment. However, for the lime being we shall take c as a constant ; application of Newton's 

law to'' the single degree of freedom 'system of fig. 'i-l then loads to the following second' 

order, linear ordinary differential equation, with constant coefficients, (ref. 2). 

- .•• • \ ...... - ; ■ 

dt dt , V 

Normally, this equation is put in more useful engineering form by inti-oducing the following 
terms: o , the critical damping associated with viscous damping (=2MaJ n >; C, the dimen- 
sionless damping ratio (^c/O; and W n the undamped natural frequency. Then eq. (M-2) 
becomes ■ . [ 

dt 

and its solution for a system with an initial displacement y(o) - y Q but^with zero j-vel- 
ocity (dy/<rt) t ^ o = 0, in an underdamped situation of C<1> is (ref. 2). 



y e 



-C<*> t 

11 [cos u>,t + — - sin UKt] € % M-M 



where u> is the "damped natural frequency" (=w n J . This solution represents a 

damped periodic function which schematically looks like Fig. M-2. 

The higher the damping constant c (or equivalent!^, C since C = c/c,) , the steeper 
the exponential envelope, and consequent ly -the vibration takes less cycles' to die down. The 
rate at which this vibration decays is of importance and, indeed, is the quantity that is , % 



1 For a more accurate model which takes the mass in of the beam into account ,, it follows from 
(\ • ref. 5, pg. 430, that the' equivalent mass in the single degree of freedom system may be taken 
" as (M + .23m). 

2 By making use of trigonometric identities, it may be showJl, that this solution can also be 

written as ' » 

j w *-\ -£w t -1 w d 



n 
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measured in this Experiment. Consider the two local maxima y and y_ ih Fig. M-2 3 in 

W> -Cu t. m <0 -TJA (i + — ) . 

which the amplitude chafes from y A a — - e to y fi * y ^ e n . Note 

d ° d 

that y B is a constant factor e~ 2irC< V w d times y . This factor is the safne for, any two 

consecutive peaks. Therefore, in general for any two adjacent peaks, v and v . the 

n ' y n-fl * 

following relationship ^olds, 

y 2ttCw /w,, 

n n d 

v " e ^ 5 

■y J n+1 

Taking the natural logarithm of both sides of this equation introduces the logarithmic 
decrement 6 . Vl-.' . . . , 

y- 2irCw 



y n+l "d 



The magnitudes of the adjacent peaks may ,be^picked off the response cyrve Fig. M-2 and the 
logarithmic decrement calculated; substitution into U-6 wi^Ll then yield £, which is a 
measure of the damping in the system. 

As will be seen from actual experimental output, the adjacent peaks are so close together 
for small damping (£«1) that it is quite difficult to obtain accurate results. A better 
technique is to us<*-4**o~peaks y^ ,and y^ +1 "^h^t^are separated by enough cycles so that an 
accurate measurement can be made. It is easy to show(ref.2), that for this case the 
logarithmic decrement 6 is r 

6 E 1 .en A" = ^ / ^ 1-7 



n y 

Finally, for small damping £«1, we may approximate eq., ('4-7) by 

1 . y l ' 



'♦-8 



C ~ Z — In 



2 Tin y . } 
* J n+1 1 

c 

A& has been mentioned previously, the damping involved in this experiment is basically 

due to the inherent internal friction property -of materials and is not associated with any 



external 'mechanisms . ' We may now write the dissipation term -c^r in eq. ('4-2) explicit 



3 ■ — ■ 

We may assume for engineering purposes that these local peaks y . and y are tangent to 
the envelope curve 

to t 
n ' n 

y = y G 

J J o to • 
d 

Comments on some possible energy loss associated with clamping procedures and also air damp- 
. iftg is discussed in the following "Experiment" section. 
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for the material damping property that wilj be measured in this experiment. From previous 
experiments on lypioal engineering materials it has been observed that the energy dissipated 
in a typical engineering material per cycle of forced . simiSoidal vibration 4s independent of 
the frequency and is proportional to the square of the amplitude. Consequently, we model the 
material damping force as independent 6f frequency and proportional to amplitude, 

• , ) , i a kr| dy 

r« matejjial damping force = - — ^- H 

where Q is the (dimenslonless ) material damping factor , and (0 Is the frequency. The 
damping force, as exprossed by eq. t't-9), is obviously independent of frequency since for . 
forced sinusoidal vibration the velocity dy/dt is directly proportional to the frequency w. 
We may now carry these concepts over to the present free vibration problem and, replace the 
damping constant c in (»»-2) in terms of an "equivalent viscous damping" based on (4-9) - 



/ 



2 J 

Md y '+ , c *L + K y = 0^ M-10a' 
,2 dt " 

dt 

MjL + M + k v - o '»-iob 

, 2 id , dt y 
dt d 



Introducing terms in a similar fashion to that used in going from (M-2) to.(M-3) and 
assuming u>^P^o> d> yields 



<*4 + 2CO) + to 2 y = o - M-lla 

2 n dt 



dt 



M-llb 



C = H/2 . 

Combining (M-llb) and- (4-8), the damping ratio C is- eliminated giving a relationship between 
H and two peak amplitudes y and y n+1 - 



n = — fcn 

mn y 

J n+1 



Equation (12) is the' relation that is used to calculate the material damping factor n • 
Experiment 

The basic experiment Involve the monitoring of the free vibrations of a cantilever 
beam carrying a heavy mass at its tip. The material damping factor n is then determined 
from eq. (M-12). It is also of interest to determine possible var/ations of n with fre- 
quency of vibration? This frequency variation is obtained experimental^ by changing the 
magnitude of the tip mass 5 or shortening the length of the beam.' It is inherent in the ^ 
discussion involving eq. ('(-9) that n is independent of frequency. (Experimental results 



5 From a study of gravity effects it may be shown that decreasing the frequency by increasing 
the tip load will yield accurate results as long as g/l < | , i-.e, the beam frequency is 
greater than the frequency associated with pendulum motion. This equation serves to put a lower 
limit on the frequency range. 



/ 

on engineering materials >at low frequency, i.e. between 2 and 100 cycles per second tend to 
confirm this; frequency independence, (ref. 3)). 

A schematic of the apparatus used to obtain data is shown in Fig. M-3. 

The cahtilever is clamped .securely to V Massive structure so that losses due to this 
connection a^e minimized. Some researchers have reported good results, (ref. ft), using a 
beam with an^enlarged end which is then clamped to a massive structure. Strain gages are 
glued to both sides of the beam near the clamped end (BLH SR-M gages wit\i Eastman 910 adhe- 
sive produces good results). The gage leads are attached to a bridge circuit so that the 
(strain produced) change in the resistance of the gage can be measured nit^ precision . Al- ' 
though self- temperature-compensating gages are the most desirable, temperature compensation 
may be attained by attaching the leads from the gages to the two adjacent arms of a four arm 
bridge', (r*ef . H) I ~ This set-up also has the advantage "of eliminate with 
longitudinal stresses (as opposed to bending stresses). The deflection y(t) at the tip of 
the beam and the strain - e(t) from the gages are both linearly proportional to tip load, 
and hence linearly proportional to each ather. The strains are then amplified and setit to 
a strip chart recorder so that a "hard copy" of strain versus time can be recorded. A very 
convenient device that combines both a bridge and an amplifier circuit, and is specifically 
built to accept strain gage input, is a Type Q Plug-In Unit for use with Tektronix oscillo- 
scopes. < This unit minimizes problems associated with calibration, balancing, linearity, etc. 
Th$ typical chart recorder that accepts the amplified signals records on heat sensitive paper. 
Chart speed is variable so that ma>£Lmum accuracy can be attained in reading relative ampli- 
tudes. Good results have been obtained using a Sanborn chart recorder, Model 322 Dual 
Channel DC Amplifier-Recorder. 

It should be noted that for metals, which typically have ^mall values of n> the 
damping associated with air resistance can be of importance. For accurate results one then 
has to consider a vacuum enclosing system around the beam. A typical device would be a long 
glass tube around the beam and base with an attachment leading to a pumj5. The vibration may 
be initiated by using a bar magnet outride the tube to attract a strip of iron attached, to 
the beam. Another "way to start the motion is by displacing the .end of the beam with a thin 
conducting wire and then sending a high current through the wire (using a battery Source); 
a high enough current will melt the wire and start the beam vibrating. 

A typical specimen is a wood (pine) cantilever 7" in length and 3/H" by 5/32" in cross- 
section. The natural frequency is about 70 cycles/second (simply let it vibrate without any 
mass at the tip and count the cycles on the chart paper). Shown in Fig. M-M are actual test 
results for- the beam with response frequencies of 16 and 30 <pycles/second , respectively (the 
chart speed is lOOmm/sec). It may be seen that it takes about 20. cycles for the amplitude 
to drop 50%; substitution into eq. ( l 4-12) leads to a value of .011 for the structural damping 
factor n- It is sometimes difficult to read the curves if the zero point is slightly off 

line, as ifl Fig* 4-**b. To alleviate this difficulty one can read the double amplitude (on both 

« 

1 

\ 

f 4 

6 * 
The actual' constant of proportionality are irrelevant since only the ratios of strain (or 

displacement) amplitudes are necessary in the calculation of the material damping factor n- 
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Fig. 4-4 Strip chart records of damped oscillations 
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vT^!«Hi4«8 of tf^jtero). Further testing with different end masses verified that this result 
V^£a rpughly ''Impendent of frequency. These results are plotted in Fig. M-5 and 'agree 
publA?^*^. Values in raf. 3. " 
T6 get results at higher frequency the specimen should be out down. Since the stiff- 
V^|y?yai3jLes a^V>>~ 3 eq. (M-l) -cutting the beam in half increases the frequency by a factor' 
'* 2 . 8 . ^he values of n for common metals such as steel and copper would fall 
i^^an order magnitude below this value for wood (thereby necessitating a partial vac- 
uumli&viroh^eyit tfo maintain accuracy) while a hard "viscoelastic" substance such as poly-. 
roetl^^e\haojfylat«»fejUcite, plexiglass) would have n about 5 to 10 times as high as the 
valsS^fet^ood . S<^t typical values of n *aken from experiments and ref. 3, are, shown in 
Fig. ^j^3)it shoula^JA pointed out that viscoelastic materials such as lucite are quite 
- -sens|&fe^>envirort%ntal conditions , in particular the -temperature. -This leads .±b._aome_. 
difficult* in reproT^n 



v -results for different specimens under different conditions, but 
it should"^ of -inter ""to the students to try to quantify these effects 
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Experiment No . 5 

Osci.l tations In a Manometer Fluid 

I ntroduct ion D ' * 

'i \ 

Ono of the problems encountered with manometers is the oscillations that can occur when 
the pressure being measured ia changed suddenly. Instead of coming to rest, the fluid oscil- 
latot, uutil the natural damping causes it to come to rest. In this experiment wo prop$s* to 
examine the relation between fluid viscosity and the damping. 

Consider the simple manometer shown 'in Tig/ 5-1. A column of liquid of length SL and 
having density p is contained in the manometer as shown. 1% the analysis will be assumed 
that both ends of the U-tube are open to air. We shall neglect the mass and damping of the 
air. For the present we shall also" neglect the damping in the liquid, ' The illustration ~shbws ' 
the fluid oscillating with the distance from the equilibrium level to the liquid level equal 

to h. • • 

The mak® bf' liquid oscillating is t & 9 

2 K X 

* . ^ ~ fc P \ ' 5-1 

whore d is the 'insicfe diameter of the tube. 

If we consider the pressure on the surface of each leg equal, the force due to the difference 
in levels, 2h , which -tcrids to restore the liquid to a level state is 

~ ' 2 r o 

2hnd pg • J " 2 



d 2 h 

Writing Newton's second law for this use, using h = — j as the rate of acceleration of the 

dt 



liquid, we obtain the expression 



T if " ' M 

or more simply 

* ^ ' h .+ 2 ? 



r\ '2 

^ ftp v h + * 2 " d P h =» ' o * 5-3 



5-M 



For more detail about the development of trhe equation and its solution see ref. 1. The 
solution of the equation indicates that the oscillation is sinusoidal with a frequency, 



n 



the natural frequency. For aj\ initial displacement h , the motion is 



h = h sin (jo t 
o n 



where u> n =. 



The liquid contributes damping doe to the viscpus action along the tube wail. Assuming *t 
this to be viscous damping, which means that the damping force is proportional- to velocity, 
h, such that the damping force equals, ch , Newton *s second law becomes 

'■ -h + - h *• |& h = 0 - J>-*> 

iji • x. - * 

For a discussion of the solution see ref. 2 or 3, and Experiments' 2 and M . Tor damping .less. 
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than critical , , the damping ratio can ba ©xpressed Jas 



o | wn p , 



whex^a c q ■ 2mw n » the critical darning factor 



(jo^ * frequency of damped systom 



\ 



(0 



88 frequency of untiampad system (the natural 



frequency) 



Referring to experiments 2 or 'I, it is possible to measure the logarithmic decrement, 
'6, and compute C"" from' i1r;--a5 - - - ■-; 



The ichoice of procedures will depend upon the amount oT damping. Tor very light damping, 

to^ and u> d are nearly the same, making precise evaluation of eq. (5r>) difficult. 

Having found C , it is possible to determine the damping factor c by evaluating 

to eq. (5-5) then c 
n . 0 

Experiment 5A * 

■ V 

Effect of Manometer Geometry on Damp ing Factor 

The objective of this experiment is to determine the effect ^ 
on the oaci.lpa.ti on of liquid in a manometer produced py 
changes in the length aiyj diameter of the liquid column. 

This experiment will involve a single liquid in an open ended U-tube. Water is suit- 
able. Now apply a pressure to one leg and release it suddenly. This will cause the column 
to vibrate. The tube will be filled with different amounts of liquid, thus changing I . 
For each & , the frequency,, co^ , will then be measured and the natural frequency 
computed from the expression developed above. From this information values of c / c c anc * 
c can be computed. Now plot, versus liquid length I , w . , c/ c c and c . Try to 
establish some i^elationsHip for c (&). 

Tor low frequency and .large damping it may be practical tx> evaluate u>^ with a stop 

watch; then to determine the damping ratio from eq. (5-7). For other cases it will be 

necessary to determine the logarithmic decrement which will require determining the amplitude 

at specific times./ Optical measurements provide'an accuracy with no disturbance to the 

dynamics. - This reJ^Jres that glass tubing be usad with. colored (or opaque fluid)- Dye can 

be added to clear fluids.. One technique is to use a photo-oleetr.i c cell, mounted so it may 

be moved vertically. When the fluid is between the light source and cell the signal will be 

interrupted. By positioning the system so that it is just interrupted by the fluid at the 

top of its oscillation, the point e^i be detected when the fluid stops intercepting the beam, 

and this height measured. Quickly relocating. the coil, and cbunting oscillations, 1 n , a 

second point, h , can be measured and 6 calculated. 
y ' n+i 

- 37 - , 
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Using a Polaroid camera (or better st;ill f . vidio tap$) wi^Jrh 3000 film, exposure aan be 
made at- two successive times. This will require advancing the film for each exposure. 
Exposures must be made ty opening the lens riear the top of the oscillation for sufficient 
time for the fluid to reach its maxinm.iv height ... ..With- a scale movlnted beside the manometer 
the height of the oscillation can be measured from each photograph. With video tjkpe the ■ \ 
recorded pictures can be replayed frame by frame to determine the maximum height at two 
•different times. Close-up lenses Qn the video camera will improve the aocyracy of the meas- 
urement. In* either case, better definition can be obtained by lighting with high frequency 
stroboscopic light. 

Another problem which ahould'be examined is the effect of fluids which wet the wall. 
Non-wetting fluids will bohave more like the model- It may be possible to treat the inner 
waXY'of tlie t ube "to "reduce 'wetTIng"'VlViri*aY9^- ~ 

An alternative experiment might be performed by repeating this experiment for tubes of 
three different diameters, and to establish the effect of the tube diameter on the several 
variables, • ■ t 

What might bo ,of more interest is to relate the damping c to the wetted area of the 
manometer tube, ird£ . Plot the results of all tests with damping as the ordinate and wetted 
s1.n7.racG area as the iibscissa. Does there appear to be a consistant relationship between them? 

Experiment SB ' 
The Effect of l iquid Viscosity on Dampin g in ^ U ~tube 

-The object is to evaluate the effect oL fluid viscosity 'on 
the damping of oscillations in a U-tube manometer. ' ' * 

Iji Experiment 5A, the objective was to determine if- damping could be' related to the 
'.tube geometry, and in particular, the wetted surface area. Our next objective' is to evaluate 
the effect of fluid Viscosity on the viscous damping, 'fo accomplish this, fluids of different 
viscosity "will be caused to oscillate in the U-tube^ and the damped frequency or logarithmic 
decrement measured. Tube diameter and length of liqirid , Si , should be kept constant for 
each of the fluids. This will establish a single value for u)^ "andwetted surface area. 
From separate measurements or tabulated values, p and the viscosity, p , can be found. 
Now determine the damping ^factor c. Ftfffcfe^ should be used with less than critical damping. 
Try mercury, and silicone oil." Select fluids for which p is known. Graphs ca£ be dratoh of 
damping factor versus viscosity. Once again, the question is:* Can a relationship between 
these properties be established? If 30, what is it? Furthermore, can a relationship be 
established for damping 1 . which includes both wetted surface and viscosity"? Measurements can 
be made with various £ and a for Several fluids in order to evaluate this relation. 

An additional exercise can be ddhe with the more visc6us fluid, That is to try, by 
extrapolation and experiment .to reach critical damping. Now compare the result to the 
expression for c , . 

»* 

c - 2ma)* 



Experiment 5C' ' . 

Modiflqatlon of t he Syst em to Increase Damping 

The objective is to introduce damping with a small 
restriction to the flow , and to evaluate Its effect on 
the oscillation of fy>id in . a U-t)ibe manometer. 

One way to increase damping in the U-tube manometer is to introduce an orifice into- the 
system. The restriction imposed on fluid movement will reduce the response of the instrument 
but Will slao provide additional damping. An orifice might^e used in the base of the U-tube 
to restrict fluid flow from one leg s to the other. In this experiment," the'restriotion will 
be placed in the opening to atmosphere on one leg. Consequently, air displaced from one leg 
by the fluid must pass through the orifice to escape or return. 

We seek, experimentally , answers to several questions. 

1. What is the effect of orifice opening on the time for the manometer .• 
to respond to a pressure change? A variable orifice or valte should be used, 
-with response time measured and plotted versus some measure of the valve 

setting or orifice opening. What ts really needed here is to minimize the 
time required to come to an equilibrium state when a change in pressure is 
made, and it is that which should be measured. It may be that too much damp- . 
ing will cause thf$ time to increase. Is there an optimum? Use water as the 
fluid. ' 

2. What :U<the effect on response time of the size of the air column of 
the restriction? 

3. What, effect is observed by changing the length of liquid in the U-tube? 

Devise a method for evaluating ' the response, in terms of these variables, and conceive 
Tnethods of .expressing the relationships among them in a clear way. 
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Experiment No. 6 ^ 

Acceleration \Qf a Spinning HUff 
A machine is to b<* designed in which a disc is to revolve in Fro* air. It is nOC69sa*y 
to determine the power required to drive the disc as a function of the speed of the disc* 
make this evaluation an experiment will be conducted uslngftb© model shown ivy Pig. 6-1. A 
circular disc in mounted on a shaft. The shaft jo supported by two bearings, with the disc 
overhanging as shown. finally the system is mounted to a firm table. 

lixperimont No; OA 

Determination of Frlctlona! Effects 

The^Vject of this experiment is to measure the frictional 
resistance on the disc and then to express it as a function 
of speed . 

The test will involve bringing the disc to some angular speed, then allowing it to 
decelerate due to the air resistance and boarTng friction. The rate of deceleration is pro- 
e, t 



portional to the torque applied by these resistive forces, according to Newton's law:^ 



Ma * 9^^) ' ' 6-1 

where I = the mass moment of inertia of the system 

about .the axis of rotation 



a s the angular deceleration 

c = tprque due to friction forces , a 

(unction of the angular speed , to. 

Inertia of the disc and shaft can be found by calculations, knowing the geometry and 
density of the disc and shaft, or it can be measured by me ants of the calibrated torsion wire 
of Experiment 1. Express it iS Ibm^ inches^ It" is necessary to measure the acccleraf i on . 
Tli.is is more difficult to do than to measure angular velocity, and it is Known that accelera- 
tion is the time rale of change o.f velocity. In terms of infinitesimals, it is the derivative 



of velocity with respect to time: 



6-2 

dt 



t 

w 1 1 e re w = angular ve 1 o city » 

t = 1 i mo , 

One soon learns that derivatives are difficult to evaluate experimentally or graphically as 
these usually entail measurements at discrete time 'intervals,,' or approximating tangents to 
graphs. It is now a quest. ion., -of making accurate measurements of angular velocity^ either at 
'discrete time intervals or cont i nuonjSly in time. 

Procedure sugg estio ns . 

Some pleasuring schemes c.ome to mind' quickly which extract excessive amounts of energy 

from the sytem. Thus, the instrument contributes to the resistance giving a false result. v 

Instruments- in this category include mechanical or electric tachometers. A simple revolution 

counter can be made by attaching a small manget to the shaft, which will pass a pick-up coil 

with each revolution ' of the shaft. As it does, an eJectric impulse is induced in the coil by 

f -% ■ 

^ the magnet which can bemused to drive an electronic impulse counter, or. to produce bissajous 

- Ml - ■ 
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FIG. 6-1 ROTATING DISC ASSEMBLY 
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patterns on an oftcilllscop© ( ref. .U and 3). 

Another simple speed-measuring technique which can bo employed is to use a Strobotac to 
bring a mark on the wheel, to an apparently stopped condition. One soon -Jjaarna it is better to 
set the frequency of tl)e strobotac and measure the time at which the wheel speed becomes* syn- 
chronous, rather than trying to measure the time-varying speed at a given instant of tijne. 
This is beat done by setting the strobotac to a frequency, say 600 flashes per minute, then 
marking the tinje that the mark on the disc appeai^s to be stationary. A single mark will 
appear at 600 rpm. There w„iil bo two marks at 300 rpin; four at \50 rpm and so on. This has 
the advantage of spacing the readings more closely in the low speed range. Repoating the 
test with the light set at MOO flashes per minute, wilT"TDTow measurement at 200 rpm, 100 rmp, 
etc. 



be tormina t ion of resistance 

It is now necessary to determine >,a .frpm the time-velocity data. This can be done by 
approximating the derivative by a finite difference approximation 

where the subscripts refer to successive measurements. As an alternative, data can be plotted 
and (angents measured. An improved technique would be to fit a second degree polynomial 

through three successive data points expressing C) (t) as 

•J 2 

(t) = A t A, t + A., t 6-M 



o J 2 



With the coefficients A^ and A,^ determined 



for a .value of t in the center of the two-step interval- J 

Now the value of c: in cq. (6-1). can be found for several values of speed w , and 
a graph c vs . w can be constructed. It might look like Fig. 6-3.. Notice that for the 

units of c (w) are in lb.f-.ft,., which requires that the inertia of the shaft #nd disc, 1, be 

• 2 2 2 

-converted to ''slugs- ft (lbf-sec -'-ft ) and a to be expressed in radians/sec .y 

Sample results 

An experiment was conducted using a solid aluminum disc 12 inches -in diameter having 
Wo 9/16 inch thick annular icings fastened with screws to the front and back faces. These 
rings had an outside diameter of 12 inches and an inside diameter of 8 inches. The disc 
assembly was fastened to the shaft of -the apparatus, rig. 6-1, and a pulley was fastened to 

sthe back end. The assembled rotating parts had a mass moment of inertia of 0.0751 lbf-ft-sec 

2 * 
or 2.M19 Ibm-ft . Using an electric motor and a v-belt drive, the s\stem was brought to a 

speed of "about* 1000 rpm and the belt slipped off the pulley. Using a strobotac, and the 

technique described, the time required to decelerate was determined. ' It is shown graphically 

in Fig ; 3-2. w - 0 . 

- - 
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FIG. 6 "2 "DECAY OF VELOCITY WITH TIME FOR A SPINNING DISC 
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From points on the curve, the acceleration wa$ computed by eq. (6-3)^ Patting I ( t 

and a in eq. (6-1) c(u>) was calculated. Results for c(u>) are shown plotted v$. spfced 

2 2 * 

IftTIg. 6-3 . For T « 0 .0751 H>f-ft-Beo and in rad/sec f is expr9?a©d in -ibfHFty- - 

■',...(*' 

; \ * ■ , * • • * 

Experiment No. 6B * . . - t. * . / 

A Iggast Squares Curve Fit" . ■ . L 

* The objeot now is to use the data fat* cAo>) from the previous 

experiment and express the Junctional relationship as a poly- 
nomial 4>f degree two using the method of least squares. * 

The form of the polynomial is 

2 

c^ ((o) » a + a^u> + a^w * S 

The method of least squares is described in refs. 2 >% 1 } and 5, and involves a seleqtiori Of 

the coefficients a , a, and a. such that the quantity, 

. ■ o 1 . i . * 

£ (a + a.'O),^ t a 0 u) 2 - c llr ) 2 ^ ^ 6 ~ 7 

k*o oik 2 k lk 1 * ^ 



is a minimum. Here the c are the values ofi c computed from the test results at the 

lk 1 • 

oi - There are a total of n such points. The method requires the solution of a set of 

K 

simultaneous algebraic equations, linear in the unknowns a^, a^ and a^ • 



s 



n n 2 



$ i n 

na + a, E u>, + a 0 * E w. = E c ,-**t/ 

,o , l k=0 k 2 k^o k ^ *=o lki 

V 

a ? w *+ . a, So), 2 t ... a ? 0) 3 ^ 8" c u) 6-^« 

o k=0 k lk=o k 2 k=o k k=o -LK * 



n 2 n 3 ' ft '* n 2 * 

a E (jQ s + a. I a), + a E u> = E c .W, 

o k=o k J k=o k 2 k=o k k=o 1* k ' . 

The procedure then is to take the computed values of and the corresponding 

o> and determine the >r best" fit according to tfte criterion of least squares. The values 

k * * > 

of ^ , a and a. y will be found using the eq. (6-8). These will 1 then be ■ substituted in 
o 1 * 

eq, (6-6) to give the second-degree polynomial having the "best 11 fit. 

Sample rgsults " " 1 : 

For the data plotted in Tig, 6-3, the expression for ^(u)) was found to be 

c^U)) = 7.6 x 10~ 2 + 7.5 x 10"" 5 u) .+ 2.05 x 10"* 0 .(0 2 

Experiment No. 6C 

Evaluation of the Effect of Air Resistance 



The objective here is fo determine tjie effect g given area 
1 change has on air resistance, and to determine hqw the 
t , coefficients of the resistance Equation (Experiment 6B) are 

affected. . 

< . • 5^ 



1. • 



Variations of Experiment 6 can bo introduced by altering* the air resistance of t^diso 
through tho addition of vanes or spoilers to the -disc , .or adding roughness to the fiat surface. 
* Suoh changes would 4 effect the air resistance which potild reveal itself by changes In o (w) . 
Another variation is Vo construct^a disc having the same m*$s and sajne moment of inertia but 
.a different diameter. In this way the bearing resistance should remain the same aa for the - 
first disc, but because the surface area has changed, the air resistance will bo changed. 1 
In this way the air resistance can be isolated from the keatd.ng friction. 

One meat\9 of changing tho disc area without: changing^fss or inertia is to construct 
an annular disc as illustrated in Fig.. 6-'« . For example, suppose the material of the annular 
dif|c has a density 2.8 times that of the first (steel/aluminum*) and t^e edge discs are 1/0 the 
thickness of the first disc, then .for each value of outer diameter there is a certain inner 
■- diameter- which will produce discs with mass and-inertia -equal to a v solid disc-, - .Values-of - the 
inner radius, r , are given for several choices of r 2 in the table .below. The radii are 
made dimensionless by .expressing them as the ratio of the radius to that of the plain disc. 



Outer and inner radii for annular disc having mass ai^d inertia of solid 
disc of radius R. Radii are given in terms of the ratio to the radius 
of the sQ^.id disc, 



inner radius 
R 



outer radius 
R 



0.9 

0.866 (3/M area -of solid disc) 

0.8 

0.773 



.795 
,723 
.blM 

0 . (^solid d isc of 

dense- material) 



Notice that if r^ is .773" times the radius of the first -disc, the second disc is solid. 
An annular disc wa£ cons t I'uct ed from steel, with an outside diameter of 10.390 inches and an 
inner diameter of 8,6 .inches. * Tor the solid disc having a diameter of 12 inches, th^js cor- 
responds to an outer radius ratio, v, ? , bi' *8&6. This results in, the area of tjic annular 
disc being 3/M that of the solid disc. 

Values of c^Coo) arc again determined experimentally and compared graphically for the 
two discs. Coefficients can then be obtained for the second degree equation using, a least 
squares criterion, then compared to those found for ^he solid disc. The question of which, 
coefficient in the polynomial is most affected is of importance to evaluate how air resis- 
tance contributes to the total torque. The change in resistance can be related empirically 
to the face area of the discs. " . * 



v 
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Experiment No. 7 

% 

Quid Friction in Plpo Flow 

■ In this experiment the objective is to reproduce the well known 
curve of friction factors for fluid flow in oircular tubes. The 

* tubes wUl he made from sections of glass tubing, or' copper tubing 
with short connections for attaching a manometer.. 



The te st facility ^ ■ ! . * 

~ 9 ✓ . . . 

A sketch- of the test set-up is shown inTlg^7-l. for thls-text, 5/16. copper" tubing ■ » 

was used having an inside diameter of 0.137 in. Pressure taps were located 60 7/1$ in. 
apart, and' there" was a section 25 Inchest 79 diameters) upstream of the first^tap to allow 
the flow to stabilize and develop. A section of tubing ^extended 10 inohes beyond the. down- 
stream tap. The taps were connected to a U-tube manometer with an entrapped air BuM>l.©; as 
shown, on whi.ch" V pressure drop was measured in centimeters of water. At the downstream end, 
the water was collected, and measured over a given time interval.- From this. the flow rate 
could be found. The tubing was clampedj$to a table and leveled. 

Water was supplied from a constant head tank. This oan be constructed easily, and 
proves to bd a useful' tool in the laboratory. A 5-gallon drutn with an : open top, such as a 
paint drum, serves for the tank. A 1-inch diameter overflow tube extends ( through the bottom. 
In the tank shown, the tube could be slid to any level by a collat arrangement. A threaded 
knurled ring was used to pull a rubber, 0-ring, seal tight around the tube; which scoured the 
tupe in place. Also on the bottom of the trfnk a tube connection was" fastened, for the water 
outlet. Water was supplied by a hose having a perforated ^leeve over the outlet to act as a 
baffle. A vertical plate with openings at the bottom acted as a second baffle to.separate 
the disturbances set up by the inlet flow from the overflow, and outlet. Final-y a gage 
g-iass was attached to permit observation of the water level in the tank. The tank was ele- 
vated on a stand to provide the elevation necessary for flow. Position of the tank oh the 
stand could be altered for head changes with fine adjustments in head made by changing the 
height of the overflow pipe, * % 

T e st procedure % ^ ' 

The level in the constant- head tar^kwas first established and the manometer reading and 
flow per unit time were measured. The level was changed to provide 'a new flow rate, and the 
new data recorded. About one minute was allowed to stablize flow before readings were ma^e . 
Flow was collected for one minute using a stop watch, and the volume was determined by 
measuring in a graduate. Water temperature was recorded in order to establish the viscosity 
of water. J 

Calculat i on procedure ^ 

Velocity, v, was calculated as 

v = % , ft/sec. 
A 

. where 0 = volume flow, ft^/sec. 



7-1 



A = cross sectional area of tube, ft^» 
- 51 - 




t 

FIG. 7-1 TEST TACIUTY FOR DETERMINATION 
OF FLUID FRICTION IN TUBES 



Ndxt~the Reynolds number was found 

«■ 

vD 7-9 
Ko * — _ 

whore D * pipe diameter, ft. 

2 * 

v » kinematic viscosity ft /sec. at the 
water temperature . 

Finally the friction factor, f , was found from the definition 

V 2g . - 

where h. « pressure drop duo to friction, ft*. 

1, * length of tube between pressure taps, ft. 
o a acceleration due to gravity, 32-2ft/sec' 
v and D ai^e as defined above. 

Friction factor waa then plotted against the Reynolds number on log-log graph paper, 
and the results compared to' the Moody diagram found in any elementary text on fluid mechaniqs 
( for example rof, l.or ■?)■ ^ 

' * i 

So me results 3 ,< 

A word should' be said about units and -dimensionality , as it is an important educational, 
aspect of this experiment. Both the Reynolds number and friction factor are dimensionless and 
their importance should be emphasised- 

Furthermore the units must be consistent, and .in this experiment some conversion has 

3 

been necessary. As Q was taken as cc/min. it is necessary to convert it to ft /sec. 
Hence, for the 0.137 in. diameter pipe, 

3 3 
r p C ,- in i f 1 ft i r l min-i 

I—. J LO-061023±Ad L 3 J L J 

Q min cc 1720 in 60 sec 

v = -~ - ~ " " 2 

T3.1M16J [(0.317) ? in 2 C— ^ 

1MM in 

_ = (0.0010739) Q, ft/sec. 

cc s 
where Q 'has the units (jjjr^O 

r. 

. _ CQ( .0010739)— D [0.317 in] — ] - 
Re = ^ • ; .. sec , 12 ft ' 

[.v <±- ] • 

* sec 



For the water temperature which ranged from 4MF to MGF V was taken as 1.58 * 10 



ft ^ /sec. 



Hence * 



Re = 1.79bMb Q [.Dimensionlesg 



r 
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FIG. 7-2 RESULTS OF EXPERIMENTAL DETERMINATION 

OF FRICTION FACTOR 
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Finally th« friction factor i» ., j ^ 



,D» f cm] [2D .^2.2-^3 - • • ^ 

"f sec * 

* ■ — — y ~ ~~ 

t L \ v ' r C0.M37S in- n r 2ft ~i r<MJ fi cm n 



- CO. 01108] .(Dimensionloss) 

\ . v ;"■ 4 

Calculations were tabulated with th<? data as 'illustrated In the,tabla below. 
Two examples* are given. 

Data k 



— j . . 



Trial 'Pressure drop Flow -Rate' - —Velocity Reynolds - Friction 

* No. (cm) (cc/min) (ft/seo) No. , • Factor 

S 2.3 596 0 .6M00 1070 0 .0622 

10 3.8 967 1.0384 1 1736* 0,0390 

A graph of friction factor versus Reynolds number was constructed .from the results. 
This is Fig. 7-2. Comparison are very good both for the laminar flow range where 

f = &i and in the turbuleht range. The data compare best to 'the 1 results for a pipe 
having a dimensiQriless roughness of — - 0.001. 



Some modifications 



Experiment 8 can be modified by adding elbows, or valves to the system and determining 
pressure dr6p as a function of the pipe Reynolds number. TJ>kse can be compared to existing 
data. • 
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EHper<A\©nt N6. 8 -* . 

Unstead y Flow of a Fluid ^ 

The object of this oxperimnt is to predict the time required 
to empty a tank of liquic} artd to evaluate the preAiotion exper- 
imentally. 

(Note; This dm be pfcrfoiAnad as an experiment without the 
analysis if it , is felt that the analysis is too difficult.) 



Introduction 



/ ■ 

Designing systems in Engineering requires application of theorie^o predict the 
behavior of the system before the system is actually constructed. It goes without saying 
tlVat "reW^ for successful designs. Furthermore y it -is inueh - 

le&s expensive to .modify the design on paper than to alter the system after its construc- 
tion. The objective of this experiment is t6 predict the behavior of a system with a 
simple theory, theh construct the system and evaluate its behavior experimentally. The 
system selected is the emptying* of a tank of water through a simple piping system. There 
are many real situations where tanks of liquid (or gctees) are emptied or filled, for 
which this experiment might serve as a model. For example, a rural water supply system, 
operating from a single standpipe might be required to supply water for* emergencies at a 
prescribed rate and for a set time. How should this standpipe be elevato4? If a hole 
were accidentiy made in the pressurized cabin of an aircraft at high alt/itude ,*how long 
would it take for the pressure to^ reach ^dangerously low levelsf? This ^nvolves flow of a 
compressible gas which is^ beyond the scope of this experiment. 

Description of apparatus 

The apparatus for the experiment is illustrated in Fig. 8-1. It consists of a 
5-gallon tank with an opening in the bottdm, to "which is* brazed a short tube, 3/8 inch 
outside diameter. Using Tygon tubing, copper tMbing aftd piping, a simple plumbing system 
is added. This can be varied but, ir» the experiment shown in Fig . 8-1, it consists of a 
90° turn in the tubing, a l/M inch (gate-type) valve, and a 1/M inch pipe elbow. 

Development of the theory . - v • 

The tank is set at a prescribed elevation above the piping. It is filled with watefr 
to some level and allowed to drain through the pluibing system. The objective of the 
analysis is to predict the ^elevation of the surface of the water h, as <| function of time 
From this information, it is simple to compute flow rate or. the velocity of discharge. 

Bernoulli's equation relates the elevation (head), velocity and pressure in the steady 
flow of an incompressible fluid without losses due, to friction. Between two stations , 1 
and 2 (designated by subscript), the expression can be written 
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where 



... . « V.T-, . 

3. -W 1 



h a elevation , ft «■ * „ " 

v « velocity, ft/sec 

% » • g * acceleration dUS gravity, 32.2 ft/sea 

p » static pressure, Ibf/ft 

Y 13 weight of fluid' pe^ oubio foot 

h '= *a constant in the flow, called "total ■h^^3f^;; r ? ^ 

Thft specific weight of the fluid is related to the density by^ Ntwto^'ai^ Sy*\** 
boliGally, , * \ ^ .» **** ; 



here 



8^ 



p = density, lbm/ft 

.._ a .*_ _dimensionless_. conversion .f actor 

° / 2 . 

32.2 Ibm-ft/lbf-sec . 



Notice that the specific 'weight is given in lbf/ft and has the same magnitude "as the 

3 • . ^ ' 

density p . For water at -room temperature, we can take p = 62. 'P lbm/ft . \ 



A -practical modification to the Bernoulli equation is to recognize "£hat flifere .. will * 
be lossed in a real flow causing the right hand side of eq. (8-1) to be deficient by an 



amount, h^, the loss is the flow. 



Then 



2 v 2 

h l + ^" + Vp l = h 2 * 2T + YP 2 



8^2 



For the emptying tank, \ake station I as the free surface of water in thft' -tank-, r -fend 
station 2 at the point of discharge as shown in Fig. 8-2. Let h ]L = h^ a the d^t^jice fronT 
the point of discharge to the surface of the water in the tank. From the same Reference , -y 
h 2 = 0. The pressures on the free surface (p^ and at the point of discharge * Xfcj)* 
arc both equal to atmospheric pressure and can be set to zero. Finally, the tank;di^\, : ., 
eter is so much larger than the diameter of the jet at discharge that the velocity "orf 
the tank surface, v^, can be* neglected when compared to the discharge velocity, v^. 
Then the modified Bernoulli eq. (8-2) can be written 



v 2 ^ 2g (h-hj 



P-3 



To bring in the time dependence of h, Write the equation for conservation of foa&{? '? 
(continuity), 'equating the flow rate from the tank, Q T> to the flq| rate from the pip^ 



V 



8-M 



For the 'tank, the flow rate can be expressed as the rate of change v of level with time,-t, 



dh A 
eft T 



ffD 2 dh^ 
C 4 dt 



8-5 



where 



A.j, = the cross sectional area of tanj^/ft 
D = the diameter of the tank, ft. 
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from the pipe, the flow rate oan b« related to t1>e velocity of discharge » Jv n , and 
* ■ ■ 

the area of the pip$, v 'A' , which with eq. (8-3). is 

-P » 

". ." ■ ' . \ % ' ' A p v 2 |t ' 4- 2 « < h -V' • ' ' °- 6 

where • " d » the inlside diameter of the pipe, ft. , 

Equations (0-5) and (8^6) can bo substitued into eq. (o-»0 to express the time depend- 
enoe <?f 1 \\ as - - « - *• 



W6 now seek an expression for h^ . For this, we will use the loss coefficients 
found i\\ elementary texts in fluid mechanics for the plumbing system. There aye: a sharps* 
edged entrance , "a "straight pipe ~ a gate ""valve I", and Harf elbow - Using "a Reference s\Soh ag 
Streeter (1) or ^Drenkert (2), these are / 



\ a t*r\ t% n t* to « 1 1 r» m • \ % 



2 ' ' - 

sharp-edged entrance; h 0 v 0-5 v^ /2g ' 8-8a 



a pipe of longth L: h^ = f (v ? 2 /y) (yd) ^ ^6b 



P 



gate valve wide ppen: = 0.-19 V ? /2fe 8~8c 



v 



90° elbow: 4 h^ * 0.9 vjy2g 8-8d 

. e " . 

The value of the* ^friction facf Or .. f in eq. (8-8b) is a f function of the Velocity of 

the flow, the viscosity and density, of the- fluid and the diameter ofth'epipe. For laminar 

flow it can be 1 written as , A . 

f . ^ .. ' 8,9 

pvd ' - " 

where ) p = density of the fluid, lbm/ft , 

v =, velocity, ft/sec . * ' 

\x = viscosity, lbm/sec-ft 

d = pipe diameter, ft. j 

Note: f is dimensionicss , and the expression pvd/|/ is the flow Reynolds number, and is 
dimensionless . A value of Reynolds number larger -than 2000 in pipe, flow indicates that the 
flow is probaBly turbulent and eq. (0-9) should not be used. -Instead-* graphical value of 
f as a function of the Reynolds number g(RG^| should be obtained ( see.. x*&f ./l/'or 2). * 
When these expressions for head losses are substituted into eq. (8), the result is' 

If . a " jf - (h " 5T + 0,19 + 9,9 f J 8 - (Re)/) ■ 8 ^ 10 

For the. case of laminae flow/ the friction factor is an easily expressible function of 
Re, eq. (8-*)), making oq. (8-10) . 

2 v 2 

* s _<L (h _-2 (1 ; 59 1^^)) ; . '< 8-ii 

dt D 2 2g • ■ d pv 2 d . , 



' • • , • • , . • • ' 

Equation (8-il) is difficult to solve because .v is a function of and appears in non-. 
— - linear form. Consequently ; th.« solution -must be a- numerical one. A trial value is foupd " 



Reynolds number. The equation is 



* * ' h ^' c 'y_ 41.59 + L /d g (Re* ) y * 



An improved e'stimatl^ofr -the - velocity -can -now be made -by -using eq. _( 8-=6) with_ the ja^ly 

computed yalue of ; , 1 ^ 



At this point, .the improved velocity v 2 (1) should be compared to v^.-. 'One way is to 



ask if 



|v 2 (1) -v 2 | *6 



where e is some small value assigned to test the accuracy of the^ iteration . If this is 
not true, the new value of v 2 , v 2 (1 \ is used in eq. (8-i2) to find a new h^, which \in ■ 
turn is used in eq. (8-6) to find an improved value of' v^ "When the iteration has con- 
verged, that is |v 2 U) - v 2 |*re, fhen eq 1 . (8-10) is solved for |£ . ■ Integration of this 
ekpressiofr^n be carried out from the initial state by . any of several methods including, a . 
simple technique called Euler's method (see ref. 3). In this. method anew value of h is 
found from the expression , . \< 

, h (1) - h (0) O (0) 6t 

/ • r dt 

In this case h (0) is the initial value of h ; (^) (0> * s the value found ^ ove S' and 
St is some Irtiall interval of time for the next step. Having found h U ,' the^ process 
can "be repeated. That., is, a new trial value of v 2 'found { h^ calculated; an improved 
value of v 2 calculated; the test made, for conversion; iteration until the best values 
are found for, v 2 and h £ ; and then Euler' s methof applied to" find h ... Gradually the 
integration goes forward until h has reached some prescribed value at which l&SRf the- , 
..computation ..is terminated.. For each step in time, the velocity and head, h,, af#known. 

It is evident tnat this is "a tedious calculation . If a digital computer is available,' 
a program can be writtdn to handle this. 'In the discussion of an actual- experiment , a 
program Written' by students performing thiis experiment is presented. ,* , v . '. f. 



, for v 2 with vh^ 'set to zero in eq!T(8-6). ' { ( „ , , "'. t 

d 2 h 

" « .. • v =.— a- 1/ 2 $ b - " • '• - "■•■'» 

• D • . 

... • ; ■ '< , . . ' • • 

(.Then a value of h. is found using equations (8-8a) through (8-8d) and an evaluation of the * - 



Experimental procedure ' . "'V-;,,-, ■ % u -f 

When the calculations have been made, a curve of head versus time can -be constructed. 

Now, using the apparatus, the tank can be filled with water and measurements of head versus A 

6 9 : . ... • %l ■ 
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time recorded &9 the tank is allowed to drain. This can be done using a stop watch and 
recording the time when the surface of the water crosses pro-selected points on a meter 
♦ stick; or, collecting the water in a graduate and measuring the time necessary to collect 
a certain amount of water , then converting this to change in heqd. For a small diameter 
piping 'system v \ and large tank diameter, the readings with the watch can be made with 
reasonably good accuracy. 

Sample results * 

In the model used by the students, the water level was rea<J in a gage glass which was 
attached to the tank as shown in Fig, 8-1- Tha. tank was filled with a clamp on the hflse. 
At time zerq, the clamp was retooved and the tiink began to drain while leadings were 
recorded. Two sots of readings wer*e made, and first' differences were evaluated to 
determine whether the' results 'were linear, (First divided difference/* are found as the 
difference. between two successive Values of head divided by the time difference between 
them.) As expected, the first differences were smaller for higher times than they were 
initially. That is to say, the change in elevation for a given time change is less for 
the smaller head because the velocity is less. This curve of the experimental data, 
Fig. 8-3, looks very close to a straight line. 

In the computer program, the students entered an average value- of friction factor, 
0.032. This simplified the calculation as it was necessary neither to calculate Reynolds 
numbers nor to use a table look-up for evaluating f . This could have introduced error in 
that the value of f lr/d for this case, with a 62 inch long section of 3/8 inch tubing 
was , 

which was the dominant term in the expression for h^ . ■* Values of -loss factors, diameters; 
and a trial value velocity were read into the program. The program, Fig, 8-M , is self- 
dxpianatoi^y , 

V| The predicted curve shown in Fig. 0-3 shows a slower rate of draining than the 

measured curve . Possible sources of t^rors arc the use of an average value for f and 
the rather large value of e used in the program. In this case, e = 1. However, the 

^Ltudents f ^approach to the problem, the qualitative "agreement between analysis and experi- 
ment , and ^h^, recognition of sources of error made this a worthwhile exercise. It pro- 
vides an opportunity^ for theory-testing and gives the student a chance to' relate analysis, 
digital computation and experiment in a single exercise. 

Modifications of the experiment * • * 

.* * 

Several modifications' of this experiment are possible. For one thing, the piping 

can be altered in. .many ways, A system of piping can be employed, and a s'ingle clement 



added or removed V ^vtand the .results coihpared . Particularly interesting is the addition of a 
globe valve to the system. Comparison can be made between the losses with a globe valve 
and a gate valve. 

Fluids can -be changed. Using a glycarine of known "viscosity will provide results 

ERIC , . • 



■i,-^t— ,~T* 




FIG. 8-5 a. 

BUCKET SHOWING 

PRINCIPAL 
DI MENSIONS 



RUBBER SEAL 
COMMPRESSED 
BY BLOCKS 



9 

ERIC 



8-5b BORDA 
MOUTHPIECE 



i 
i 

n 



, 1 1 



i 



8-5c SQUARE 
ENTRANCE 



65 - 




8-5 d BELL- 
MOUTHED 
ENTRANCE 

vr' 

ih 
i; • 



± 



r 



8-5 e SECTION OF 
BELL-MOUTHED 
ENTRANCE n 



<5 x 



?3 



T 



r r 



which are quiCa different from experiments with water. 

The height of the tank can bo changed m can the diameter of the piping to insure 
turbulent fiow. 

Using a pressurised tank equipped with a pressure 0 gage , the time for a gas to 1 leave 
the tank can be measured. Here t the important term Jle the tank pressure expressed as a 
function of time. This experiment should include" an assumption of compressible flow. 

Any of these can be done as an experiment* without introduction of the analysis 
which may be too advanced for some first and second-year students. 

Experiment 8A ' * - ) 

The Effect of Entrance Geometry - - - * - - - 

in this experiment ^eyerai entrance conditions will be employed 
and head losses determined for the draining of liquid from a 
tank . . " ° 



Introduction * - 

The previous experiment was modified, to investigate the effect' on the emptying rate 
of the tank for several inlet geometries. Three were used. Eaoh is illustrated in Fig. 8-5 
The first is a Borda mouthpiece,. Fig. 8-5b. In the second , the v tube has been lowered to 
forjn a sharp-edge entrance Fig. 8-5c. The third was a bell-mouthed entr^ice, Fig. 8-5d.. 

A five-gallon drum was used for the tank with a 'l/M in. diameter brass tubing. The 
tubing entered through a rubber seal fastened to the base of thfe tank. This is shown in 
Fig. 0-5a. . The tuning could be pushed through the rubber seal so that the tubing end was 
set at any selected distance above the bottom, or'flush with the bottom. A smdll bell- 
mouthed entrance was mad£ of aluminum. It had a 7/8 inc)} diameter Opening at the entrance. 
The bottom was recessed to fit snugly over the end of the brass tubing as shown in 
Fig. 8- be. 

Referring back to eq. (8-2) and the analysis that follows, it is apparent that thex 
only difference between the previous case and the present one is that h ♦ is not zero, 
but is equal to the length of the tube. Consequently eq . (8-7) becomes 



Sr =- -r V2g«M. 2 > - V 



L 2 

D 5 

As the tube length is fixecl, h is a constant for all of the experiments.' 

j» 

The experiment 

The tank was filled to a pre-de termined level with water. ^A Stale was fastened 
Inside the tank to provide a measure of the level of the water surface while the iowex^ end 
of the tube was closed with a stopper. One of the inlet conditions was established. 
Several minutes were allowed to pass for the fluid to come to a quiescent state. At a 
signal, the stopper was removed and a stop watch started. At a succession of closely 
spaced intervals, the time was recorded and, simultaneously, the scale was road to indicate 
the ll&ight of the water surface. 




FIG. 8-6 HEAP -TIME CURVES FOR A DRAINING TANK SHOWING THE EFFECT 

OF THE VORTEX FOR TWO EN TRANCES. 
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FIG. 8-7 HEAD LOSS FOR FOUR. INLET CONDITIONS 
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Plotting h V8. time, the- values o£ ■ — Qould bo determined for a set of valuer of 
h. This can be done graphically or by approximating the derivative in finite difference 
fofm 



dh^ h (2) -h (1 > 
dt ^ — 



Here h (l) and h (2 ^ are values of h at two times, separated by a small interval 6t . 

With d, and D - measured from the experiment, and jj- were known; a calcu- 

lation for h. was made fl^om eq. (0-13). * 
Graphs can bo prepared for the several cases of , ( . h vs. t and for h fc vs. h, to 
compare the effect of each entrance on losses When the experiment is run, it will be noted 
that a single vortex is formed which is quite spectacular. It becomes apparent that its 
presence affects the results, and better evaluation of the entrance conditions, can be made 
if the vortex is destroyed . This is possible if several oieces of window screening- are 
put in the tank dividing the tank into pie-shaped sections. The screening "should go from 
the tank bottom to above the water surface. A modification. of the experiment , which is 
id only with this vortex, is presented as Experiment • 8-B. 



concernec 



Sample results . 

Measurements were' made during emptying lor a tank having a diameter, D, of 11.13 

inches, and a tube having an inside diameters ,. d , of 0.55 inch^Four entrance conditions 
were used : . ■ 

Condition 1. Borda entrance, 6 inches above base ' 
Condition 2. Borda en franco, 3 inches above base • ( v 

Condition 3. Bell-mouthed entrance, 3 inches above base 

Condition '». Sharp-edged 'entrance — tube moved to a position flush .with the 
bottom of the tank . 

. Tests were made without screening (vortex present) and with screens (vortex absent). 
Graphs of the head, h, versus time were constructed for each case . Shown in fig- 0-6 
arc conditions 7 and 3 with and* without screens'. Trom these data fi- was approximated 
.and the head loss determined. Head loss,- h £ , is shown for conditions 1, 2, 3 and M with 
screens in Pig. 0-7- 

These results indicate that the emptying times ar«,not appreciably reduced by the 
elimination of the vortex. Also indicated is the conclusion that the losses are measurably 
less for the rounded entrance than for any of the other entrance conditions.' There is .not 
a large difference .omong cases 1, 2 and M . Erratic behavior in the values of h £ is due 
to difficulty in evaluating the derivative accurately from the loss-time curves. 

It is probably more useful to plot, the head loss in"a "dimensionless form, one which 
compensated for the change in velocity ..wit* head. . Using an expression like the friction- , 
factor or f has this characteristic , * h » 
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FIG. 8-80 SECTION SHOWING WATER 
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FIG. 9-8b SECTION SHOWING 
DIE TRACES IN THE VPRTEX 
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Her© the Velocity can be defined as an .ideal velocity, if losses were fcero, , ., < 

* • " ' 

or.^es the measured velocity, - 



'8 -is 



^ „ Q - » f^'l) 0-16 



/ - v 2 ^ ~~T "* ~2 ( dt ) 

<nd , d 



examination of ©q. (8-16) shows a sgirce of error in t^iis experiment that might be' 
excessive. When . the pipe diameter, d, i small compared to the tank diameter, D, the 
head is easily measured iti the transient Is- it is changing slowly However, the ratio 
D 2 /d 2 is a very large number and must b* multiplied by a .small number, dh/dt, ^ which is. 
difficult to obtain. It is always difficult to determine derivatives from d^ta. for a 
curve of a given function, the curve of the derivative is less -smooth than the curve of the 
function, and the integral produces a curve that is smoother than that of the function. 
Crrors in data are magnified in differentiation and suppressed in integration. Conse- 
quently, for this test, more accurate results might be obtained by measuring 1 Q directly, 
this can be done by collecting flow in a graduate for a short .period;, finding q as the 
amount collected por second; and computing v ? from the expression 



v 



? " A 2 

■nd 



This will be the average velocity over the interval taken to measure tho flow. 

equations can then.be loft in terms of v Q and the losses computed as before. 

Experiment 80 

/ * 

Study of a Vort ex ( * 

The objective is to examine a vortex in a tank being 
drained. 

As indicated i\ the previous experiment, when flow was established in the draining 
tank, a vortex was formed. The profile of the water surface is shown in Fig. 8-0a. An 
experiment can be devise\to study the hehavior.of this vortex under different flow 
conditions, and to trace tHe fla^tTt^' vortex system. 

The tank is filled wife- water, and allowed to settle for ten minutes. Flow is begun 
by removing .thVplug at the -hot ^ of the p.i, P e . It may be found that the vortex will not 
; -devero P ;..s'oon enough ifVthe water iWoo still, and the tank could empty before the vortex 
is established/',. A .small swirl can vbe Wouraged before "draining to induce the vortex. 

Food" coloring cah be adjected through a long hypodermic tube at different points in 
the flow/l!o\.ob^erve%t'!re direction of f low \c*fif ferent levels. Surface movement can be 
deleted bV'«*v6rtex. meters'"' ; . These cambe Constructed from small bits of paper, with a 
crOsVoh .^ph,; which' Tloat. on the surface. . 



Sample results ' . 

The.vortox meters will show the suriace velocity in the vortex and indicate whether 
the surface flow is "irrotational" or not. It is found that thp velocity decreases as one 
mpyes away from the canter of the vortex . The motors rotate loss in the region^ away from 
the wall and not: loo clos& to the center^ If therein no rotation of the motors, the flow 
^Ts behaving as an ideal free vortex. Viscous effects cause rotation In the Vortex which 
will' cause the meters to revolve on thoix* bwft ayes. 

Observation of the dye traces shows the entire mass to be swirling. The dye will be 
/-seen to move in along the bottom of the t$nk V and then spiral upwards alor^g the outiet pipe 
as shown in Fig. 8-8b. When it reaches the inlet, it is di*awn Into the tube, 

i " References 

L, Stree.ter, V. L. , Tluld Mechanics, McGraw-Hill ^ New York, New York (.1971),; p. 9a7. 

2. BronKert, Elementary Theoretical _^_l u i ( i Mechancis , John Wiley and Sons \ New York, 
New York (i960), p. 22M . " ~~ ' 

3. Arden, I}, and K. N. Asti.ll, Numerical .Algorithms , Addison-Wcsiey , Reading, Mass. 
(.1970), p. 221. 





Experiment No. 9 
Convection Heat Transfer 



- When "an objeot tft.one temperature is immersaa In -a fluid .. at another * temperature . 
heat'is transferred from or to the object. The mode of heat transfer Is called conation. 
One measure of the rate of oonvective heat trans fer* is the coefficient o^heAt' transfer ^ , 
h. It Is defined from Newton's law of o.ooling as * • 

w ^ erQ • q- = rate of heat transfer per unit of time ^ 

1 . .- -A ' k. 3urface .aro.a- Q.L tlie Qbject { 

t = surface temperature of the objecj 

t = bulk temperature of the fluid,. ( 
- ^ o i ' 

*» « i 

•lYpical units of h are Btu/hr-ft 2 -°F. Sometimes h is grouped with a length dimen- 
sion, ,W and the conductivity of the fluid medium, k, to .form a dimensions number 

called the. Nuaselt number, , - 4 

hL „ 9-2 

If the object is in a gravitational field, and the fluid movesVer ft>e to gravi- 
tational effects, the heat transfer is described as natural or free convection. For the ^ 
same fluid, the rate of heat transfer can be increased by blowing or pumping the fluid past ^ 
the object. This is forced convection. One soon learns that the magnitude of h depends 
on the fluid velocity, the geometry of>he object and /luid properties. «ioh as conductivity, 
•k; viscosity, \i density, p ; and specific heat, c p . 

in this series of experiments wo shall determine the coefficient of heat transfer 
between a small object and a fluid for several situations by means of a transient temperatur 
measurement. The objects will be small enough with high conductivity' (copper or .aluminum) 
so that we can assume the temperature of the object to be uniform, bet that temperature 
of the -object be t. From eq. (9-1), the heat transfer between the fluid and the surround- 

in8S 13 - q = -h A (t-t Q ) 

Jhis will be reflected as a change in temperature of the object. For a mass , m, with 
specific heat, c , the rate of heat transfer is delated to the time rate of change of 

P • • 

temperature by the equation 

. * dt 9-3 

a - m c -f— 

• ■ P dx 

where T is time. . 

He can now. /cord the temperature history of the object as it is heated or cooled. At any 

point in time we can obtain t and yjr • 

. / , ' * o mvi' a few the obiect and then solve fo^the coefficient 
We can obtain values of m, c and A lot tnc oujcoi . « 

P vBBfj 

of. heat transfer by equating the tvo expressions for q, W • 

.dt, 
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Experiment OA , « . 

Determine the Coefficien t of Heat' Trans fer for Natural Convection 
r ' ^ . ' . 

• The object is to determine the coefficient of heat transfer in 
natural convection from a email copper object by making a 
.transient measurement of the temperature. 

' The object to be tested .is -a small cylindrical copper slug with a rounded- end. This 
is hollowed out and' mounted on a plastic support ' as sTu>wn. inTig. . 9-1. Some dimensions are 
give* of>the figure. Embedded near the' center of the object^ is a copper constant!* thermo- 
couple with the loads carried out. through .the plastic rod- For information on thermocouple, 
see repSl or 2~' A diagram of th* test Tet-xrp is given Tn Tig.;9-2: (For -natural- con vectio 
eliminate" the air supply and hot-wire anemometer. ) • 

The copper slug .is first heatec\.by putting it in a blast of hot air, or immersing it 
' 'in boiling watep. If: water is use.d, dry the slug before:' beginning the test. Using the 

t" laStic support and i ring kand, the slug is placed in still air with its axis horizontal, 
he temperature 'is then measured at specific 'intervals . A better approach is to use a. 
millivdlt recorder such as a Leeds and Northrup Speed^ax recorder.' . The chart record gives 
a'curVe of^mrUivolt. output of the thermocouple versus time which can >e easily converted, 
to- temperature v f us time. A chart record is shown in Fig. 9-3. .Air temperature is read 
from a second .thermocouple suspended in -the free- air. ! - 

'* Derivatives of temperature with -time .can 'be brained from the chart record and the 
coefficient of heat transfer computed. at various temperatures along\he-cooling curve. 



' ' V tor the- copper probe' il lustrated with c > 6,092 Btu/lbm-°F., the product. 

* > 



mc ^ © 

— P- -~- o.ii3M ntu/ft - r 



A • o 

At a. temperature- of 'l77.3F, the change in' temperature was measured as 5°F in 37 seconds N 

from 'the cMart. From -this;, ■ " j .; ■ . 

. . • "' • .' . dJ_ ^ _ %- (-a600)- u =: -i|86.5°F/hr. 

- dT . V ■ . ' « . 

With the air temperature measured- as 78F,'the coefficient of heat transfer from eq. (9-M)'. 
. ' i * ' . ' ■ 

. . • u - (O.U3M)(U 86.5i . 2 12 -Btu/hr-ft 2 -«°r. 

, *" " • (17'/. 5 - 78) ' 

Converting 'this to a Mussed number usin/as the char^t eristic length the dieter of ' 

the copper.slug, 0.7S in/and k = 0 .01-b But/hV-tt-V -for air , ' s '? ' 

. . ' ■ - Hu . = (2.12) (JS* 8 . 8lt * • ' - . ^ 

Curves of/.usselt number versus V temperature can be developed and^ thefts compared 
to those .given in ref. 3 or M . . »" - - 
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Modification* 



, In natilrAl convection- the geometry is an important factor. The experiment* own be mod- 
ified by making the measurements with the axis * vertical," and the "blunt noge pointing down- ^ 
ward. A % second Slug made of copper but with a diameter of 3/2 inch was used to evaluate- the' 
effect of diameter on the convective heat transfer. Results for \v were about 10 pbrcwU 
higher for the small probe which is about - what w6uld be ejected for free convection from 
horizontal cylinders (ref. 3). To evaluate the effect* of , properties *of the slug on the 
experimental result's, "an aluminum slug carr'be made having the same dimensions as the large 
probe. Differences between boat transfer for an aluminum- and for tlje 'topper slug vkzvq less 

than the observed errors in any given test run. . .. v. . , . 

Instead of air, the response- in 1 an ice bath or hot water can bo measured. This 
requires a recorder ^To obt ain the temperature history, as the temporaturci response is too „ 
i*apid to follow by manually balancing at potent iometer. ^yS' 

Experiment . 9B • ,j 

Determination c^the Coefficient of Heat: Transfer in Korce^ Convection 

- Xhe object* is t<* determine the heat transfer coefficient , • 

in forced convection by makingpfcransient measurement- of 
/ the temperature in a copper object. 

This experiment. differs from tho farmer In that -the oL^ect is placed in a stream of 
moving air. ' * Air* can Ve supplied by a commercial vacuum cleaner, a woman's hair drier or a 

.gun. The object can bo heated/or .cooled and put in a stream of unhealed air, or it 
cdh beat rqpm temperature and placed in a stream of heated air. In any case, it will b<3 
necessary to measure the temperature of the air being blown past thto object, using a ..thermo- 
couple suspended in the free air stream. It is also Important to be ^\ble to change the 
approach, velocity of the air and to measure its speed just ahead pf the object.* This 
should be done with the object in place, bu^rfot necessarily during the time that the 
temperature record Is made. A\imple hc^-wire apparatus capable of indicating the magni- 
tude of "the air velocity works ^jesCT^ 

The slug is brought to some temperature different frctoi that of the air stream. It is 
then placed in the stream and the thermocouple regord made on a x strip chart recorder. Using 
these results tcs> evaluate t and dt/dT, eq. (^i)'is solved for h, and .the Nuaselt num- 
ber calculated. Glyphs of h versus slug temperature and . h versus flow velocity can be 
made*. 'Most useful are graphs of N.usselt number versus the \^tow Reynolds number- Here the 
Reynolds number is defined as „ \ v 



Re = 



pvd 
M 

where *■ < ■ p = density of the »ir, lbm/ft 



X /( p =s viscosity of the air, lbm/sec-ft 

v - velocity of the air, ft /sec 
d = diameter of the copperas lug , ft 



'I)imensionl$$s correlations can then \>e made # lining tho for 

. ■ • . - • , * sr 
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Fig. 9-3 Strip Chart 'Record for Temperature Hi^fbry 




wh«r« A and. w arft constants found f^om the plot, 
Sample roaulta 

Th« test set-up ia/tfhown in JMg. ^^r-^Air was supplied from a vacuum olenn&r 
exh&udt, with the speed of the motor controlled ..with"* VariacV The spfcfcd of the blower 
wfjus sot and the air Velocity end temperature measured ^t the -point where the slug was * to 

be placed-^ The sltig was/ heated in boiliftg water, withdrawn, dried and s4t in pla<^e jdth 

* / ) ' * 

a ring-stand support./ J t wa? positioned with its axis parallel to the air stream as shown 

in Fig. . 9-?. The time temperature record was made as done in Experiment 6A. 

Results were computed from the records 'and graphs made of h- ' vs. tefyperatui-e .and 

h vs. velocity for velocities ranging from 1000 ft/min'to 2S00 ft /mi A- Nusselt numbers/ 

were, plotted against Reynolds numbers on log-log graph v paper for several temperature^ } as 

shown in Fig. 9-M . %t is apparent that, thp effect of velocity is more important, than thp 

effect of slug temperature, on the heat transfer. . Tests were run f'br the large and small 

copper (3/'» inch diameter) (1/? inch diameter) slugs, and the results correlated for a 

range of Reyrvblds numbers from 3000 to lObOO; as .follows:' 

f Tor the small probe, Nu - 0.1MM Ro°"^ G \ / ^ ■- 

:\isa - 0/669 



For thra large probes, Nu * 6.151 Re 



Mod i f ieatiorvs 



One obvious change is to place the probe so its axis is normal to J he air velocity. 
Another modification is to induce . turbulence ahead of the object by placing a .wire cross- 
wise to the flow. In a scries of experiments using the large copper slug*, an incx^ease of 
five to ten percent was measured in the heat transfer coefficient with the tripping wire 
creating turbulence. 

Both experiments 9A and 9B offer an opportunity to consider and j whenever possible, 
evaluate sources of oxpowjmcntal error. * VP 

Referents 

l. t Schenck, I!., Theories of E ngineering Experimentation , Mfcraw-llill , New York, Nc*w York 
( 1968), pp. 9^-160. ~~ . ' ' . 

2.. Beckwith, T. G. and N . L. Buck., Mechan ical Measurements , Addison-Wesley , Reading, 
; Mass.- (1969) , pp. i»63-«»70. . . : 

3. McAdams, W., Heat Transmission , 3rd Ed., McGraw-Hill, New York, New .York (19b l »), 
pp. 176-177. r . . . 

14 . Rohsenow, W. M* fiind II, Y. Choi, Heat > Mass and Momentum Tran sfer , # Prent ice Hall, 
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Experiment No . t .10 
The Magnus Efrf«ot* 



■ . A"; . **' 



? The ob'joctive is to determine the Effect of 
■ rotation on the trajectory of cr cylindrical body 
moving through air. 



* \ r 



ft'.. 



It is well known that the trajectory of a spinning ball or cylindrical rotor 1? olir^d;. ., 
A carve thrown by e baseball pitcher is one example- The spin of the object sets up pressure *\ 
forces normal to the flight of the object' These foroos cause the path to cur?/e. This is 
called the Magnus effect. , The Fletlner rotor ship was designed to utilize thi*^ priiripipla ' 
as a mean$ of propulsion (ref. 1). A simple experiment can be performed to illustrate this j m 
phenomenon and evaluate i\s influence. ^ 

Experimental model m f 

The model was a simple cardboard tube such as used for a pap«r towel or^oilet paper . ". 

roll. In this case, the roll 'was l\ inches in diameter and 'I inches long. /\ circular* disc 
.was eoroented to "each end. These were cardboard, 3 inches in diameter and 0.025 inch Wick. 

They wore used to mark the flight. The assembly is shown. in Tig. 10-1. ' J 

.Propulsion <of the tube was acc6mplishod by using nine rubber bands tied" together. 
.These were attached to a string 50 inches long. The tube was wrapped with the string and 

piace^.orizQntally on the table- The string was connected to the rubber 'bands whiob'were 

stretched as shown in fig,. .10-2. When the tube was released, the bands pulled the string 

forward giving the tube spin and forward velocity, A little practice U necessary to 'control 

the flight. <. ■ 

Experimental pr ocedure - . 

It was necessary to heebie to measure the position of the object as -a function of time, 
to determine* the trajectory, velocity and rate of spin. A Strobotac might be used. The 
technique used in this experiment- was to ignite a corner of one of the circular discs before 
starting tl\e iligftt By 1 photographing t:he-'#a4ght on a single exposure, the trajectory as 
well as- the spin cap bo recorded . Photographs made in^thisway are shown in Fig. 10-3. 
These were made on Polaroid 3000 film using a Speed Graphic camera with a ieris opening of 
M.S.- The shutter' was held open from 1 second* before flight until 'after the flight was com- 
pleted'. Photographs were made in cts^rkened room With the camera 15 feet from the flight- 
path.' To 'record time, a vibrating reed was . used , formed by\'astening a cardboard flag on 
the end^of a steel bar which passed in front of the lens. Shown in fig*. 10-»» , this caused 
the dark regions on the photograph of Ih^Pajectory. The flag was placed two inches from 
the lens , 'and vibrated at the rate of 3.8 times per second. *A scale In the flight plane was 

' photographed to give a scale of length tcAhe photographs. It is important . that the flight . 

• be confined to a plane parallel to the. film surface.. This/is'; accomplished bettor in 
'.Pig. 10- 3b th#i In 10- 3a, 




"Suggested by ^.f esfeor ( Lloyd 'BrSfethen, Dept. of Mechanical Engineering , Tufts University,, 
Medford, Masse^jnus^ti?. ... t \ 
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FIG 10-2 VI EW SHOWING POSITION of COMPONENTS IN EXPERIMEN TAL SET-UP. 
'— " — JUST BEFORE RELEASE OF TUBE - 
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Fig. 10-3. Phdtographs of flight path of Spinning 
• ' . _■ Tube 
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Analysis of dqta and- sam ple results 

From the photographs, ,the velooity J angle of trajectory and rate of spin, oould £e 
measured.- This'was done for each of the .time segments', which were of 5/10 seconds duration. 
It was noted that while velocity was Blowed during flight, the rate of, spin did. not Quango 
appreciably . The ^variables are the spin rate , N; • th e> velocity, along' the ' path , V; -the 
distance a.long the, path,' 3; the a ngular velocity of spin , to ;> the angular direction, of 
the trajectory, 9; .and ^he time, T. Several interesting graphs -and crosarplotp can be 
made for combinations of these- Of most WeUst i^ the lift ooef ficient , (y This . is 
defined as the ratio of the force nonndl to the trajectory, (the lift) L, divided by y the 
4 force equivalent "to the dynamic pressure, ,pv 2 /2g ? . , Using the projected the cylin- 

der to relate pressure and force, , 1 "*" " 

v »■ 

r s _ ~h— , * ' io-i 

■■ 'where - = tube length , ft 

\ r *' tube radius, ft • „ * 

' . 3 ' 

p a air density , lbm/f t . .. 

. " v = velocity determined from photographs,, ft/sec - 

> • 2. 

a » conversion factor, 32.2 ibm-ft/lbf-sec • 
°o - - — . 

' It is necessary' to compute the lif t ' force:, from, the .photographs . , Consider the .sketc* 
" ".Fig, 10-5. 'The lift force is normal 'to the trajectqry . Acting on. the bod> are|.he. forces 
..causing centrifugal and tranguntial acceleration and" the .weigh t force , ..(<fig7£ 0 ) . The tangen- 
tial force .is drag and does [not contribute to the lift. -Hence, lift is equal t<y the "c*rx- 
trifugal force" te^pjgrt; th* pompotfent qf weight in the direction of lift; or . 

V. \ L s 2V> !!!& cos 0. . ' ■ - ' w-2 

*o R 8 o . ' , 

v ... . " : ' • ■ 

wher( i ^ n\ = the mass of the roll, Xbm 

6 V defined in Fig. 10-5 ' 

g * \the acceleration duo to gravity, ft /sec 
& *= j^ss conversion factor, 32,2 lbm-f t /lfrf-sec 

" ' V V ' * R >= the radius of curvature defind^ in l&g- jtO-5, ft. 

The angle 0 can be measured' from the photograph using drafting machine or equivalent.. 
To find the radius of curvature, recall.. * ' ' - • 

or . * R ~ ' . / . ,-. 

The approximation to the derivative J|' can be fquncT by taking tangents to -a agve of! S ' . 
•versus'. G, of using- a finite difference approximation" 1 from measurements o» ^ibmA Q.- 
•tlith .ft,. 6 and v determined Trofa the photographs a*^ described , and the massif* the qaper 
roller determined, lift 6an be computed from eq. ( 10,-2). " Measuring A and^ r, .and 
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FIG. 10-5 DIAGRAM OF FOR CES ON TUBE .» 




obtaining the proper value for p, the lift coefficient can bfe found from e<j. (10-1). 

Relation can be investigated between lift *cc^ffi^itefft and a tangential/ velocity 
Reynolds number or \ Spin velocity Reynolds number by graphing the results; Since th* lift 
is related to the circulation \ it can be demonstrated that|a relationship exists betttean the 
<Lift and the rdtio of spin velocity, (O r, to tangential velocity, • v. . Hyra u) is ther" 
spin rate of thfe tube, radians per second, and . r p is the tube radius. This result is given 
by Prandtl, ref. 2. , 

Tl\e lift coefficient was plotted against ■ U)r/v and compared to the -result in reference^ 
2 for tubes with discs' on the ends . These' curves appear in Fig* 1Q-6V Thp results of the 
experiment are ir> "(rery good agreement with the curves from Prandtl* 



References 



T 




Fluid Mechanics, McGraw-HiO.1 , New York, New York (jLWh) , p. '»26. 

2.' Prandtl, L. and 0. G. Tietjen's, Applied Hydtfo and Aeromechanics , Dover Publications, 
New York, Nev^ York (193H), pp. 8'4-85. ' , , .'•'.„ 



1. Streeter, V. L. , 
A 



^ 




SPINNlNGxTUBE WITH, END PLATES^ 

,, r —4 — 




S 



Experiment Hq» 11 

. * t* ; , . w The Arohimedea Sorgw Pump *' 

* The^ objective is to investigate the prinoiple of » \ , 

the' ArQhitoedes Sorew and to .evaluate ita perform- _M±_-iL 
anoe as a pump. '«••«. 

Archimedes (287-212 B.C.) is usually credited with inventing the water screw pump 
which bears his name (ref. 1).. A relatively simple model can be- constructed by wrapping, 
transparent, flexible plastic tubing in a spiral around a cylinder as shown in Fig. 11-1.* 
The tubing is open at both ends, with one end in a tank of water. Hater is carried up the 
tubing as the cylinder is rotated. In each coil, the water is trapped at the J>ottom, and ^ > 
moves up the spiral because gravitation keeps it from being carried over the loop, 

^Performance var iables / * - 

— - ; . - ■ \ ^ ^ . , 

It is evident that, if the angle of the shaft, t», is made, too large ,■ depending upon 
the" angle of the -spiral, 0, -the water in the spiral would run out and the flow rate would 
be zero., .This would occur when one," leg of the spiral became horizontal, or O = 90 - 

Design variables include: angle of elevation, a; helix angle 0; tubing diameter; 
diameter 6f the helix (shaft diameter j ; and length of spiral. ^Operating variables a|e: . 
depth .of liquid in the., tank, shaft spied,' properties of liquid such as viscosity. Experi- 
ments cart be developed around any of these Variables . As- one normally builds a single 
system, the tube diameter, shaft diamete'r. and, length , are fixed, but it is simple to change 

elevation angle or. helix angle. 

" performance parameters include the. yolumetrio pumping rate , Q; ^the volumetric 

efficiency, n„ 1 ' : and the pumping eff iciency n p . Pumping rate is easily measured by 
collecting liquid over a given time span.' Volumetric efficiency is the r 9 tio of volume ' 
actually pumped, divided, by a maximum displacement . Maximum displacement differs depending _ 
upon the velocity oV the pump. * If the velocity were very slow, the tube would fill to the 
level of the water v in the tank; that is, everything ' in the angle* 2w-0 illustrated in j 
Fig. 11-2. Tbis is called static flow volume. When the tube is rotated , the maximum yolume 
is equal to the . volume "swept out by the mouth of the tube. If the last turn of the helix 
were at 90° to thf axis of rotation, these two would be the same. As it usually isn't the •• 
swept-out volume Q is related to 6 and the helj* angle, 3, as 




Q_ 



= 'a [(R + r)(2TT-e)J to Cos3 U ~ 1 



max o 

where * r = tube radius 

R = shaft radius 

A = ' open area of tube f irr , 

o 

( (d^ i angular velocity' of shaft 

'fi and 9 are defined above. 
Then the volumetric efficiency is * ^ 

'* n = Q/Q " ' - x ^ 2 

. v max 

The pumping efficiency, n . is the rati^ of .the work done to raise the water to the 
work required to drive the screw, or 
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where 



2 . 



p, density of Liquid, lbm/ft 
g ■ fctfayitationai acceleration , ft/sec 
g « conversion factor, 32*2 lbm-ft/lbf-sec' 



Q = flow^rate, ft /min 

h* » total lift or vertical distance from water 
surface to shaft centerline at point of 
delivery, ft' * 



W = 



In order to compute .^p* 



work., Ibf-ft/min. + y 

■i 

it is necessary to measure the worlc In addition to Q. 



One method would require converting the motor to a -dynamometer . To do this, the motor 



is supported in 'bearings so the casing is free to revolve, as shown in Fig. 11-3 (see 
also ref. 2). It is constrained by a spring scale fastened to an arm on the casing. 
.The torque required to drive the motor pulley is. the product of the net force on the 
spring scale, F, and the" arm length, 1. For the scale to indicate net force, it is 
necessary to counterbalance the static torque due to the weight of the arm. The net work, 
W, is then ! , - 



W = 2TTJINF' 



ll-<+ 



where N is the motor speed in rpm. . With units bf ft in feet and F " in lbf, units of 

W are ft^lbf/min as required in eq. (11-3). Losses due to the drive system will be 

included in W measured this way. Other methods of obtaining W might include using a^ 

transmission dynamometer (ref. 3). '"This requires using a pulley with "elastic beams" for 

'{he spokes. Using strain gages fastened to the spokes, the output' of the strain gages 

can be calibrated against torque. Measuring pulley speed, the work can be found using 

eq. (11-U). *A crude estimate of\work might be obtained by driving the system through / 

lowering weights. Speed control is lost- this way ahd. results may be inaccurate unless 

kinetic energy of weights and- system is accounted for. 

* ** 

The experiment % ' ' 

A test set-up was built like that shown in Fig*. 11-1. The shaft was <*2 inches long' 
and 3}i inches in diameter with \ inch diameter plastic tubing. Ball and socket joints 
were used to support the bearings. The bottom -bearing was able' to carry' the thrust load, 
because of the flexible joints, the angle 'a could be changed. Drive was through a 
V~belt, using, an A.C. motor controlled wi$h a Variac. Water was used and collected to 
determine flow rate. Water was supplied to the tank at a rate to keep the level constant.' 
Tests were run changing variables of shaft angle, a; -depth of water in the, supply tank; 
and shaft? speed.' Analyses were 'developed to determine volumetric efficiency. and pumping 
efficiency. Numerous interesting questions can be developed concerning the operation of 
this device. Does it work if 'the end'of'the tubing is completely submerged? How dqes 
viscosity affect performance? What causes the work losses?. , ^ 
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FIG. 11-2 SECTION THROUGH SPIRAL SHOWING ANGLE 
SWEPT OUT BY OPEN END OF TUBING 




FIG U -3 SIMPLE DYNc/mQMETER 
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Sam ple results ( 

* In Fig. 11-H, a graph Is shc^m of the flQj* rate versus speed for a number of operating 
depths in the supply tank. Tests wgre made at various speeds for a fixed tank depth. The 
water. level was changed and the tests repeated. Results can be graphed for H v ahd ' 
against speed /tank depth or inclination,, ou From the results given here, the students con- 
cluded that for an inclination of. a= '*0°, optimum shaft speed was 2U5 rpm with enough water 
to cover one-half of the first coil of the helix. 
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Experiment No. 12 

A Solar Water Pump 

x . \\ .■ . ; . .• 

The obieotive is to design »->touild and evaluate the per- 
formance of ,i model solar -water pump. . 4 J. "~ "'. ' .f 7 ' 

Operating principle ' .. . 

In re^ons where power sources are not available but There is much sunlitfht, one aeeks 
to utilize solar energy to accomplish useful tasks. One proposal is a water pump using 
solar energy. Jta principle is to use the expansion and contraction of a confined volume 
6f air to pump water,- energy for the process coming *rom solar ..radiation. 



The system is showmia Fig. 12-1, A sealed volume is connected to a water source, and 
a water delivery line. is trapped between the water surface and a sealed cover plate. 

The plate" is the radiation receiver. It must have a surface which absorbs solar radiation 
efficiently" and it must conduct this energy into the ai* bubble.. Heating or cooling the air 
will cause the pressure to rise or fall- * ; - 

The operation can be described by a cycle. Solar radiation raises the air temperature 
at constant volume. This is shown as process 1-2 on the pressure volume diagram for the air 
(Fig". 12-2). At some pressure,: p 2 , a one-way check valve in the delivery 1# is forced 
open and the air expands, forcing the water through the delivery line to the storage tank. 
A second check valve in the suction line prevents water flowing" back through that line. As 
the delivery head, h is essentially constant, this process,' 2-3 in the cycle, is ideally 
at con stant pressure, p^ At point 3, the air is cooled. This could done naturally by 
waiting until night, or shading the surface. It could^be done more frequently if S/>me of the 
water from the tank were allowed to flow cvef the- plate. At any rate, the pressure will drop 
during cooling, and with both valves cldsed, the process, will be at constant volume I - 

When the pressure, , is reached, the suction-line valve will be opened anM water will be 

f drawn into the pump chamber. As suction head, 1^, ' is constant, the process, 4-1,' will 
, ideally be at constant pressure. Atmospheric pressure would normally lie between ?1 and 

p as shown in Fig. 13-2. 

v There are several performance quantities of interest. In the actual device, the vol- 
umetric flow rate, Q / the total head developed, 1^ + h 2 , apd the working efficiency n p . 
are important: A thermodynamic cycle analysis can be developed^ evaluate the thermodynamic 
efficiency for an ideal machine. This we call ^ and it can be shown to be a function^ of 
the pressure ratio, p^ an 4 the temperature ratio, T^. These temperature are the < 
maximum and minimum" attained in the cycle and are expressed as absolute .temperatures . 

> ; , > 

* Cycle analysis „ 

To develop the cycle efficiency , we require the following concepts: 
1. Assume the trapped air obeys the gas law-* 

— * =* constant # 12 " 1 

where p, V and T are pressure, volume and * 
' temperature respectively. 

- 95 - ' " . 
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FIG. 12-2 PRESSURE - VOLUME DIAGRAM OF SOLAR PUMP 

CYCLE 



2 . Wopk <Jone by a gas during a aonatant pressure prorata* 
between states ^nd B is 



i' 



'3'. Work done by a, gas during the constant volume prOce^i 
•" - la xero. •/ - 4 ' ■ . ■ 

H. Heat transferred to the gas in.a obastant pswiuf**. change 

- ^"from A to B* is . : ? y ■ 



Where c is- specific heat at eonstant pressure > 



5. peat transfJpped in a oonptaat vplume change ; *H i«S ~ 

. - ^ * H (T« -T s ) - - 

v v B A « ; 

> \- . * 

t Where c is r specific heat at constant volume\ 



6. ' For a, thermbdynamic c/cle, the suimnation .of^the w<^^0 
equals the summation of thereat transfer. 



7. The cycle efficiency is the ratio of the ^ w^k^dpn^ 
divided by the he$t transferred, into t»i^as , Hj£ 



or 



Work 



12-2 




r . ,12-3 

'I ■ "V • 

1- ^ 

12-4 



12-5 



As the net work can be expressed as the heat t^ansfei^d 
into the gas minus the he^t transfer f romv t^*, gas , H^, 
the efficiency can .be written as ^ f 



H 

o 

4q 



The net work of the cycle can be expressed with eq. (12-2). Defining this quantity 
e's* W. ' the jiet work >*is done -along paths 2-3 and 4-1 in Fig. - £2-2. Hence, 

) ' 



12-6 -ft* 



noting 



W. 



V - V 

V 3- 2 



p 2 <v 3 ^ v 2 > + Pj (V x - y, 



. In. terms of the'pump head v 



W, 



V - V 
. 1 



(p 2 - p x ) 



g o ' 

- a h 



6V» this becomes* 



■'V 



it*** & 



12-rl 



where 



1 "»j 



p = atmospheric pressure, lbf/ft 



2 



density of liquid pumped> ibm/ft 
acceleration due to gravity, ft/ 
g, = 'mass conversion factor, 32.2 lbm 



&nd .h. and h are suction and pumping head, ft. 

v 1 * -2 ii . 




12-8 



■ ' ". -To express efficiency we can evaluate H and H . Heat is transferred .into the 
... cycle 'along^'ths 1-2 (constant volume) and 2-3 (eonstant pressure). 
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Th%n 



H, 



°p < T 2 *l> 4 % (T / V 



12-9 s 



Heat is rejected (transferred from the air) Wng patb« 3-»4 and iT-1. The absolute value 



of H is 

o 



/ - 



V (T 3 .V Vp/V V/ 



Putting 1?heee exp^w^iotis' in the efficiency eq.. (12-6) 

After some algebraic manipulation and expressing C p /c y * k, this equation becomes 

. ■ V T ' " T T *' 



1 '- 



T T 
3 2 

T~ — — T — 

!3 



-T, l;r ~ 1^ k ^ (1 - ST-)] 
1 l l 1 3 



12-\0 



12-lA 



ltf-12 



Now we 'can define two cycle parameters which are Important . These are' the pressure ratio, 




R 



p 2 /p l 



12-13 



and the temperature ratio, 

T R = 

In addition, noting that 
equation (12-12) becomes 



n i 



T 3 /T l 



T T 

T ' T 
1 2 



R * 



= 1 



T P T 

, /[P R -l + k/(f 1)] 
*R K -R R 

1 1 + k p r 4: 1)] 

1 *s « 



1 



12-15 



12-16 



* Pump capacity is important and related to the volume change 6V. Applying, the gas 
law, eq. (12-1)' to the'constant pressure change, 2-3, we can write . • 



V 2' 



5V 

1 + — 
1 



12-17 



Using the expression, eq. (12-15) for lyTj, we can write 

T. 



6V 
V, 



R 

R 



\ 



12-18 



An expression for efficiency can be developed in .terms of 6V/V JL and P R by substituting 
eq. (12-18) into eq. (12-16), So we have expressed the important performance parameters for 
the thermodynamic cycle, n. and flow (SV/V^ in terms of Jwp dimensionless ratios T R 
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A aeries of calculations was made to evaluate the cycle efficiency in terms of Tj^ 

and' P using eq. (12-10). Calculations were made" for three temperature ratios, 1.1, 
1^25 and 1.50. It Is doubtful that temperature ratios larger than 1.5 are attainable, - 
Curves are shown in -Fig., 12-3 . It is apparent that the efficiency la zero for.V zero pros- • 
sure ratio r Efficiency is also, zero when the pressure ratio equals the temperature ratio. , N 
At this time, the temperature at .2 has reached the value at 3 so no ; oonstant pressure process 
ia possible. The oyole is. completed , by returning to point 1 along path» 2-1 with .ho work 
being accomplished. Consequently, there is a -pressure ratio .which gives a well, defined max- 
imum for cyole efficiency tor any given ratio of temperatures . It would be well to design 
for that . condition . . 

The -'tes t model ......... _._ .... . _. ....... . * 

A model system can be built ^ and tested. One is shown' in Fig. 12-4. Itf .Consisted of 
a steel tank 18 inohes in diameter having a capacity of 3 gallons, A delivery tube and ^ 
suction tube were brazed into the bottom. These were 1/M inch diameter oopWer. - One-way 
check valves were installed in each line, The top plate was pf 1/8, inch thick aluminum, 
held in place with 34 bolts. A rubber ring gasket was used to seal the plafte and tank. 
Gasket cement was applied to insure a good seal. The tank was checked for /leaks by blowing 
air into it while it was submerged in water.- The suction pipe was put in k one-gallon 
jug of water and the system was primed. , 

As a heat source, two 375-watt "heat" lamps were used, placed' a few 
top surface of the plate. The plate surface was" given a dull black coating. Thermocouples 
were located in the air space and in the water below. 

•Cycling was accomplished by turning the lights off for cooling. Thi 
adjusted so that the system pumped water, and the cycling begun. The twcj 
sibie a maximum temperature of over 200F for the air, which was reached ijn 35. minutes . 
Although a higher temperature might have been reached, the temperature -had begun to level 
off,- and this point was selected as point 3 (lights off). This can be seen in Fig. 12^5. 

• ' ' * * ' ' i J'.' 

The test ' » r # ♦ 

. ■ i 

It required at least one cycle, to have repeatable results. During. each cycle, the 
water was collected and measured at the delivery, and the amount -taken in during suction was 
noted. Temperatures of the air were recorded every minute, and for. the water every five 
minutes. Cycle time was . 70 .minutes , with 35 minutes for heat ir^. Times were marked at 
valve openings and closings ., Results for tv*o cycles are .show^in Fig. 12-$. Here tempera- 
ture histories are given for the air and water. Along the bottom are noted the events taking 
place. "H'« stands for heating and "C" for cooling. Numbers correspond »to states in the 
ideal cycle. > "• 

\ 

Sample results 

For the test described above, some measurements wejr>e: . \ 

tank volume * 3 gal. or 11,353 ml. 

" water volume . ! *216 ml. \ \ t 

air volume 7137 ml. ; \\ 



inches ;from the 



•delivery heVid Was 
lights made pos- 



,.V|. • 
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roo» ai*; tati^tfatura 

/ Vol uma p umpad par cycla 

(avaraga) 

»* 

Head, h 0 



Head, 



AtTMpph«Hc pressure 



lQ.in. 

29.88 in Hgf 



From these data we can determine pressures (p=62.4 l^ri/ft ) from eq. (12^8)T 



32 2 0 5 

p 2 2-9.85 (.491) + 62.H (^y) ^ . 



/ 



• « 14.66 + 0.22 « 14.88 psia. 

^ 32.2 , 10 



Therefore 



P x 1*4.66 62.4 
14,88 . n _ 



144) 



s 14.30 psia 



P' = 
R 



■ -a 



From the te^t ( Fig. 12-5), T 3 207 + 460 * 667 R and. 



73 + 460, = 533R 

V SSI a 1 - 251 

, R . FIT 

From th<* curve (Fig. 12-3), the ideal cycle would be expected to 'have an efficiency of 1.2%. 
In this' test, the heat flux to the plate was not measured, so no comparison could be made. 
Using eq. ( 12-18) » the flow per cycle for the ideal case can be determined as 



1 = 0.202. 



In the /actual cy^v'ihe volume of water pumped (which equals the volume change of air) was 
^OHO^ml. For an original aii^Vo Lyme of 7137 ml., 



>^ v Test 



10«*fl 



a HO. 116 



7137 

This is son*wfiat less than the p#a*i,cted value. 



Work and efficiency 

For the model, or prototype ^mp ,M* *P>?* M?? ful "° rk » W » can be ex P rrfS8ed in 
•terms of the flow per cycle, and the head /or ' . . / 



W 



Q (h + hj, ft-ibf/cycle 

SL 1 2 



12-19 



For a measured initial air volume, V A , Q is related to the volumetric ratio as 

5V ' "• 



Q ■ 5V 



The fenergy supplied to the jlevice per cycle, E, must be measured from the incident 
^r^diation on the surfsce> For a lamp source this' might be difficult to get. For a solar 
source, tabular values are "available for ' radiation from the sun normal to the sun's rays 
for. various latitudes and seasons. An average value of 425 Btu/hr can be used. Knowing 
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the angle of inclinatioh 1f)WtHpeh a normal to the plate and' the sun (a) the radiation 
received per cycle is ' _ 

• H ^ "(Area of plate) (Solar radiation) cos a (6t) * 12-20 




' where fit is the hour fraction for hbating. . ^ 

• While it is doubtfigti|hat solar Radiation on the plate could be as 'intense as that 
from the lamp, let's make jjjji estimate of pump efficiency assuming that the sun's radiation 
j£"Hpplie<j the energy in th<* [M^t . described . Pump efficiency will be defined as 



' 12-21 



1 

The volufme pumped per Cycle '<kks 



> tips' 

From eq, (12-19) 



Q =" *0'«0 mil/28316 mJt/ft 3 = 0;037 ft 3 



°%; ( 62 '32;2 ' '. 2 ) ( ' 5 + ' 833) 

3.0dMt-lfb/cycle or 
3 V .08 Mj r lbf 

778.(ft%f/Btu) = °' 0 °39,Btu/cycle 



Assuming the sun's inclinatipn as M5 .V^nci noting the heating cycle took 35 mitmtes, then 

v x \- - v 



from eq. (12-20) . , v.. \ 

* i .v ■ 

■ A *2 * jsfe 

. , H R = U25 (cos ^5)(^ 1 i y^ ) (||) ' % 

' \ ■ ~ ' 

= 310 Btu/cycle. , ' « 

* . \ ' . ^ 

Finally, from eq. (12-21), the actual efficien9y is. 

p 310 

This is considerably less than predicted by th& thermodynamic ^analysis . There are, losses 
due to heat transfer to the environment from the plate and from the air to the water. While 
the temperature rise of The water*was small, it \nust be remembered that the heat capacity 
is very large compared to the air. A flexible* partition with thermal insulating properties 
could be used to separate the air and wate* thereby reducing this loss/ In spite of the 
low ideal efficiency and .the even lowei? actual Efficiency , these are not. critical when you 
consider that the energy is Vfree" . 

Modifications * 

~ . < . ,* ; 

Tests can be made to evaluate design variables such as: 

V 

a. changing surface characteristics of the plate, 
, b. varying the initial air/water volume, • 

c. varying the delivery head h^ or the suction head \\ • 

d. evaluating the effect of. water cooling on the surface, 

e. testing under actual conditions (solar radiations) . 
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Experiment Np.* 13, 
mtxaigi4g-^^i- natlon of M ^«t>ial Properties 

prepare* by Robert Greff 
Associate Prpfeseer df tteotfanicai Engineering 
Tufta University s 

. V.' f 



Tha object of thin axperimant is to_uae ultrasonic wavaa to 
datarmina proparties of engineering materials In an efficient 
and accurate manner. „ 



Introduction 



In thia work, some elementary concepts in wave propagation and elasticity theory ara^ 

■ • ^- * • * * 

combined so that common 'material properties yean be found experimentally. Tests are. performed 

» by transmitting shqrt bursts of high freque/oy Vibrations (ultrasonic wajte trains) into mater- 
ial specimens and making calculations based on the measured "times between reflections. Since 
the specimen length is, known, these time intervals may be easily converted "to wave 'speed; the 
relations between wave speed and material properties are th#ft used to find results. Typical 

r * 

material properties determined in this manner by students ax^e Young r s modulus, shear modulus, 
. and Poisson's ratio and their respective variation with temperature. ^ . 



4 



Nomenclature 



E 

G 
V 

P 
I 
t 



c . 

o 



Young's modulus s 0 r 

shear modulus (= E/2(l+V) for isotropic materials) 

v 

Poisson f s ratio 

m r 

mass, density - 1 * 

length of specimen ^ • 

round trip transit time for extensional wave 

round trip transit time for torsional wave 

velocity of extensional waves 



Basic principles * 1 

* / 

The basic principles of ultrasonics involve using an ultrasonic transducer to convert 
electrical energy to acoustical or vibrational energy (ref. 1-3). In the experiments to be 
described,. the reverse process is also used, i.e. conversion of vibrational energy to elec- 
trical energy. Ultrasonic waves are usually defined as high frequency waves above the limit 
of human hearing, i.e. above 16,000 HeTrtz (cycles per seafcnc}). The frequency range generally 
used impractical ultrasonic equipment ranges from 20 kHfc "(kiloh^rtz) to 25 MHz (megahertz) . 
Two typical- methods used to generate ultrasonic waves involve 'the piezoelectric effect and 
the magnetostrictive effect. -» 

Piezoelectric transducers m^ke use of the fact that a varying electric field can expand 
and contract certain crystalline materials. Some typical materials that exhibit this peizo- 
electric effect are quartz, barium titanate , lead zirconate titanate, lead metaniobate, and 
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lithium sulphate. Applying* vojLt^ge across a piezoelectric material will lead to a thick - 

S * 

ne3S variation proportional to the voltage variation. Conversely, applying a force of a 
certain frequency to this piece-y^ll generate a voltage of the same frequency. 

Magnetostrictive transducer^ make use of the* fact" that a varying magnetic field can 
expand and contract certain ferromagnetic and f errimagnetic materials. Some typical materials 
that exhibit this magnetostrictive effect are nickel, remendur, and permendur. This effect 
can be reversed; a mechanical stress applied to a magnetostrictive material causes a change 
in intensity of magnetization." ^ 

In a typical ultrasonic ^application utilizing the pulse echo technique, waves of a 
oertain frequenpy are u^uallyjproduced in either of two J^ays ; in the first method, the fre- 
quency^of the u applied electrical current is set at t&e desired l^vel while in the s^cohd 



method a shock pulse is used, with the frequency then determined by the resonant characferis- 

tics "of the transducer. Typically, the transducer is driven for a short period of time 

result ^g in a short burst of ultrasonic waves. The pulse travels * through the material to M 

the^opposite boundary and is reflected back to the source as an echo (Fig. 13-1).. After the 

transducer, has transmitted the short burst of waves, ?t then acts as a receiver for the 

returning echoes. The echo signals are then amplified and displayed on an oscilloscope so 

that the time interval between two successive echoes can be measured accurately.* 

v 

Experiment " , , 

The object of the experiment is to obtain material properties for engineering materials 
using low' cost ultrasonic accessories. A typical experiment that will satisfy these require- 
ments involves the pulse echo method for thin w ire s pe cimens (ref. as shown in Fig. (13-2). 

The (iltrasonib equipment that a typical laboratory would require are .a pulse generator and 
an ultrasonic transducer; the total cost of these items can be purchased for less than $600. 
For example, this equipment can be purchased from Panametrics (yfaltham, Mass.) as follows: 
$395 for the pulse generator ( "Panapulsar" ) and $195; for the ultrasonic transducer ("Modulus 
Transducer 11 ). General purpose pulse generators are also available from manufacturers Such as 
Hewlett-Packard (model 21UB, approx. $800) or from General Radio (model 1217C pulse generator 
and model 1397A power amplifier, total cost approx. $950). All of th^e pulse . generators are 

i 4 ' 

used £y just plugging into a* regular 60 cycle AC wall outlet . without the need for any addi- 
tional power supply. 

As shown in Fig. (13-2), the pulse generator energizes a magnetostrictive transducer, 
launching ultrasonic waves down the lead-in ^wire^. The specimen** is joined to the lead-in 
wire by using bonding, brazing or welding. The waves are reflected from the front and from 
the back of the specimen. These reflected waves (echoes) travel through the lead-in wire to 
the. transducer which now acts as a receiver. Because of the magnetostrictive effect these s 
waves generate a voltage which »is shown on the oscilloscope. The time interval, t , between 

Enough time must be allowed between transmitted pulses to permit the return o.^ the desired 
echo ^signals aYid the decay of any subsequent reverberations. '» 

The Specimen should be thin compared to "Wavelength ; the largest cross-sectional dimension 
should be less than 5mm for a wavelength of 25mm. f 
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FIG. 13-1 SCHEMATIC FOR PULSE -ECHO METHOD 
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L. 
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L. 



ENVIRONMENT 



FIG. 13-2 SCHEMATIC FOR THIN-WIRE EXPERIMENT 
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these ybltagfc signals is aqual.'to twice the v time & for the wAves to travel through the sped- 
men.^ (Typical pulses may be rGad off a calibrated oscilloscope to an accuracy approaching 
0.1 psec). T^is timd interval leads to the d^erminati&n of the modulus/ 

To determine the modulus , we may use some results from elasticity theory (ref. 5). For 
an ed&stic solid which ffi thin compared to the wavelength of the propagating waves , it may* be 
proven that Young's modulus, E, is i 

E = pC 2 \ * 13-1 

^ o 



where p is the mass density ai\d C q is the velocity of propagation of e^tensioiwl .^waves . 

If the* length of the specimen is then the time interval t is V 

J - ° , * 

2£ 

t = ~ . 13-2 

.. _ - . o_ .4 C . . ._ 

( ° . 



so .that from (13-*) 



4p (A/t ) 2 13-3 

o 



Equation (13-3) is used to find numerical values of modulus. It should be noted that the 
length *of the specimen is relatively unimportant. * Lengths from 10mm to 10 meters can be 
employed with a typical length being 50mm. 

■ Shown in Fig. 13-3 is an actual photograph of the oscilloscope trace for a test of 
work hardened stainless steel 30Uwire. The, diameter of the lead-in wire (which is also 
SS30U) is 1.5mm while the specimen diameter is 0.75mm. The lead-in length is 33" (8Hcm) ' 
and the specimen Vength is 2" (5»,08cm). In Fig. 3 we see a double exposure in which the 
top picture has a (horizontal) time scale of 50ys/cm while th^. lower trace has been expanded 
to lOys/cm and is therefore much easier to read. The echo pattern shows the time between the 
-interface echo (i.e. the interface between the lead-in and the specimen at which there is an 
abrupt impedance change associated with the diameter change) and the end echo (i.e. the 
echo from the end of the specimen). This time is t = 19ys, and from eq. (13-2) tt>e wave 

4 O 

velocity C is 

o • , - 

.* » 

* C 2£/t ^ 

o o 

* = 10.16/19 13-U 

= 6.53 cm/ysec 

6 

-which may.be written as C = .53 x 10 cm/sec. Then from eq. (13-1) it follows that 

2 

E = p C 

° ■ 6 9 

~ • = (7.8) ( .53 x 10 ) 13-5 

12 ' 2 
= 2.2 x 10 dyne/cm 

= 32 x 10 6 lb/in 2 > * 

T 

In this manner, material property data can be obtained for a great variety of materials 
including metals, ceramics and plastics. Some, typical results obtained by using this 
pulse echo method are listed below. x 
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INTERFACE ECHO— 1 r-END ECHO 
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INTERFACE ECHO J Lc 



END ECHO 



FIG. 13-3 PULSE ECHO TEST . FOR STAINLESS STEEL 304 
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Table of Young'a Moduli (lb/in x 10 ) 

i Aluminum r * 10.4 

Bkr^liuny/% \ ' r , ■ - N • ^ 

^ Copjfcr -"^ "' 17.0 <f 

" r Nickel ; > 30.0 ' * > .j " 

J Lucite 0 .H8 / 

' % • 

* There, are aeveral-jnodi float ions that can be made In thV preceding experiment £o obtain 

additional property d^ta. Sbear modulus and Poieaon^-ratidv v may be obtained by using 

a. transducer that generates .botlv extensional and torsional Waves simultaneously (reff 6). 

_Since torsional wave speed i&^about >*& of ex^ensioneil ifftvfLjspeedi foy many: mater lalg. t the tim# T 



separation between the two types of waves may be controlled by varying the length of the lead-- 
in wire. It is possible y by proper impedance matching bfet*(een the leAd-in wire ancj the speci- 
men % to have the extensional,. waVe echoes damp out before the torsional wave echoes appear on 
the scope. The resulting trace on the scope will ^en .show both these wave patterns clearly 
and the appropriate time intervals can be picked of ^accurately . The shear modulus G mAy 



ybe determined in terms of the time interval 1 t between shear echo peaks using ^ 

s / 



G = 'HpU/t ) 



2. 



13-6 



Pois^on's ratio may be calculated from the relationship between E and 6 for- ah isotropic 
material - •• 



G ~ 



Using eqs. (13-3) and (13-7) V 



2(l+v) 

may be solved in terms of t 



13-7 



and 



t 5 « 

s ' , 

v = — f - l 

2t- 



13-s8a 



13- 8b 



It is interesting to note that the specimen length i does not appear explicitly in this 
relationship for V. 

The variation with temperature of the foregctfng material properties may also be 
determined by modifying the basic* experiment . The additional experimental apparatus needed 
to^include thermal effects is sketched in Fig. 13-2 in the fotfm of gashed Tines enclosing 
the specimen. Essentially one cap build a smal^oven facility to enclose the specimen and 
sub ject- it to various environmental -conditions ^The experiment is^run in' the same. f assort 
as before with time differences between peaks being read, off the sgope*; conjsetjueiitly v the 
variation' of modulus with temperature can be determined ^ It /should be tfotfedthat in most 
.solids the sound velocity decreases as the temperature ^increases which leads to ^n. increas- 
ing value o£ transit time. * - " ■ J . * , < 
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;} Enclosing part of the lead-in wire has no effect on the transit tfcme, since only .time_ 



differences are' found. 
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; Ao knowlaatteiftant • 

*' •fr4^h<^' : i>oknoiiUdg|« aaveraJL. h«lpf«>_. ^•^••ion with L/ C. LynnwortH, Department 
•V'^ Hfaii^d^trUPti^ Test inland Measurement, Panametrics , Inc., Waltham, Mass. 
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Exparimant Mo. 
Flow of Air Through g Nojgla 



* 



WW 



Tha objagtiva is to avaluata th* performance of, a nozzle in 
compraaaibia flow by maaauring tha reaction force in a 
noezla. 

K 

Note: * Experiments 1« and 14A Involve compressible fluid 
mechanics, whioh may be too advanced for students bslow • 
the junior or senior level. 



We have a small converging nozzle approximately 0.1 inch in diameter which we intend 
> to use for a small air turbine. It will be used to sxpand air from a high pressure to a 
v ' pressure near that of the atmosphere in order to produqe a high velocity jet. Turbine 
> blades will convert the momentum in the jets to work in a small turbine wheel. - Tor the 
* 'turbine to be efficient we want an efficient nozzle. -To aid in- the design V« would like to. 
i,* r < *SJhave some quantitative expression of the Efficiency of the nozzle. ^ 
jr-- ■■ ". a nozzle velocity coefficient will. be determined as the ratio between the measured , ; 

velocity of the jet to the ideal velocity. . , - 



.ti: 



s f:v..;-. .. 



V 

T 
IV 



m-i 



V^Ci 4 "- v'. ■ » . 

^ By ideal velocity we mean the velocity which should have been attained by an ideal gas 

.V flowing £^pm the high pressure to the v exhaust back pressure in an idealized way; i.e. 

^^^^^^re^si\>l)^^h6^t heat transfer. This is also isentropic, or at constant eptropy. If we 

•V" 7- ' coni^^^^fW of air through a control volume consisting of a large diameter Chamber 



) j-1 and tji&ftoj&la, we can write, for steady flow, 

fif 



2 -2 
v - v 
2 1 



Q <* ^ x K ' h 4. 

• A 2 ~ 1 7 2 



r". 



+ g(z 2 - z^ 



m-2 



w = 
h = 



v 



•It" 



heat transfer- 
shaft work, 
enthalpy , u + t?V^ 



,^or^e p^eVer^t^aSe, there is, i 



v =" velocity (hence ^— is kinetic energy) 

gz <= potential energy; product of elevation and 
acceleration due to gravity 

iti = mass flow into or out of control volume 

u £ b internal energy . 

pV * product of pressure times specific volume, 
or flow wbrk 



s 



no shaft work done; heat transfer is small enough to 



ne^leJV; change potential energy is zero; and the entering kinetic energy — i s 
sra^ll ^ompared to .(tha^t leaving . This last assumption is validated by the large 
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area ratio at the two orkss auctions. Consequerftl^y 

Here r and J are conversion factors. To convert Btu 19 ft-lbf , J is used, 

o . ■> • » 

J = 770 f g^k- ; and to provide proper units for mass we introduce * 32.2 1 ^f, aeo ^ ■ 

TheiV with h in Btu/lbm, the units are correct • In terms of units the rigtit hand^ side 

^ . 1 ' i 1 / l»n-ft . ft-lbf . Btu. 1/2 J , ft* 1/2 ^ ft_ . 

of the velocity equation Is ( y — *. ( Jr see whic " 

^ Ibf-sec ^ sec 

are proper units 'for veloci^. For an ideal gas *( assuming constant specific heat) 

h = c^. With this the velocity can be expressed in term6 of temperatures 



v 2 - .v^jc ^^ : t 2 ) i/2 

' * 223- 8 Yc (T^ - T 2 ) 



1M-M 



This does not require that the process be ideal. It is not practical to "measure in 
the high speed flow, so we seek another way of determining the actual velocity. However, 
if the ideal T 2 is available, then with this in\he velocity equation wa can find iv^. 
In terms of inlet states and exit pressure, for an ideal adi$b$tic (no heat transfer)- 
expansion of an ideal gas having constant specific heat, ^ 

it ~ T (— ) k 

■ ^ P l • 

where k is the ratio of specific heat c /c . It is simple to measure the temperature 

* p v * , . 

of the gas at the low velocity of state 1. From these equations, and the measured - 
p aYid p^ it is possible to compute iv^ under the assumption that p 2 = barametric 
pressure of surroundings. 

The objective of this experiment is to -devise & means of measuring v^ and compare 
it to iv^ as a function of pressure ratio P^/P2* 

Experim ent al s cheine_ » \ ' 

One approach to the determination of the actual value of v^ is to use Newton's 
second law, and Jby measuring reaction of the nozzle, (or impulsive force x>f the jet) cal~ 
culate v 2 - A device is to be built like that of Fig. \LM-l. y The nozzle will be^ mounted , 
into' the end of a large (1-1/2 in. nominal dia'. ) pipe. In fadt it is simple to form a 

a 

converging nozzle ip a p ( ipe gap. Air will be supplied to the large .pipe with a very 
flexible rubber tube. The pipe is then suspended by wires to form a long pendulum. Use 
wires 8-10 ft-, in length for this.. When air flows 'through the nozzle, a reaction force 
Will be set up and the, pendulum will swing to a stable" position, causing the wi_re to, 
traverse a measurable angle 0 . The force -of the reaction, F j will be related to the 
weight of the pipe and attachments, w. Using a force balance thjJP?elation is 

F„ = w tan6. . . 1,,_6 

K » . . / 
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FIG. 14 - I NOZZLE SYSTEM 
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FIG. 14-2 NOZZLE ASA CONTROL * VOLUME 
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Using the momentum' equation, applied to the control volume (the pipe) the fcipce can,, 



be related to velocity. \ 
In the x-direotion the momentum equation for steady flow is 

E Forced on control volume ■ E ^momentum flux out ~ E momentum flux in.^ 
With + assigned to forces'' at the right* we see from Fig. m~2. \ 



E Forcos ■ Ftp (A p - A) - P aurroundlng A - - A (p^ • 

Because the x-component o^l the entering velocity is zero, the momentum flux entering 



is zero, while the momentum flux out of the control , volume is 



t, 2 

* V 2 * ^2^ V 2 ( Qee n ^te 1 below) . 1H-7 
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where the mass has been found using continuity . A good approximation to p is — - 

• ' • • - p 2 l/k ' , ' • ' 

p « p (-£) • 1 , iM-e 

p l w 

Although this assumes reversible flow, the tfrror^ in calculating the density will be very 

small for a nozzle. Then the force of reaction on 1>he walls of the control volume, which 

is equal and opposite to F> is 2 ' 

P 2 Av 2 * ■ v 

* ~ , F R A (p surr ~ P l> " ~ i V 

& o , 

(Notice that the force is' to the left in Fig. 14-1 as F_ will have a negative value). 
We could write in terms of absolute values 

A 2 
P 2 Av 2 < 



I F„ I ^ |A (p, C p ) + , I (see note T>e*oW 1H-9 
' R' ' K l x surr g * 

o 

By measuring | F R I> P x > P surr > T i and A » we can now « c^P 1 !* 6 v 2* F R is measured b Y 
determining the swing of the pendulum and the weight of the pipe. A pressure gage and 
thermocouple in th4' pipe allow p and T 1 " to be measured. Ayea A is simply measured 



Note : 



li When the ratio of the nozzle exit pressure' to the nozzle pressure,' 'P surr /Pi 'I s i ess 
han a certain critica.1 value, the f]/tf* is restricted by the properties at ?He nozzle throat, 
hese are called critical properties, which can be designated p*,* p*, T* , etc. For ajLr, 



than 

These are called critical properties 

pVp, = t 0.53* Consequently the reaction force must be written 

c 

|r R |-|A(p i; P*) + f-v 2 | 

and m = p ft ' A v*. , 

J- • • - \ 

But v* - kg RT* 

o •* , 

hence for air, ■ - k = l.M, R. = 53.8 lbf-ft/lbm- °R , 

• • Pi A 

in = 0.532' -~ lbm/sec. • / 

, 1 . . / 

for ideal flow to the nozzle exi^. For a more complete discussion of this see, for ; 
example ," Shapiro, A, "The Dynamics and Thermodynamics of Compressible Flow" , vol. 1, 
The Ronald.-Press, New York, N. Y., 1953, p. 85 
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as 



is environmental pressure . P 8U ^ r . Now both V ? . and ly, 2 oan bs found along with k^. 
This oan be Repeated for various nozzle pressures, p^, and a graph of k n vs. iv 2 can 
' be- prepared. Pressure p. 1 oan be controUed by ^ 

5 i N 

Sa mple rggultg \ % * ' 

f ■ 

An experiment was performed using the, tost facility illustrated in Figr* 1*1-1. ' Air 

was .available from a tank suppled by' a Small compressor, The air system was capable of 



sustaining a pressure of 35 psig at *th* nozzle . 

Results for two tests are given in the following table':. 



r 



1 Deflection of 
pendulum, 9, 



degrees 



5 
6 



T °F 
75 



Nozzle Conditions , 
Pressuro' 



Temperature 



P^ P«ig 

30.0 : 
U0.O 



Atmospheric pressure 1H. 7 psi 

Weight of nozzle, piRS. an d attachments = 8.5 lbm 
Diameter of noaale at * exit •0.125 in 

Area of Nozzle at exit\ = n( >125) = .0.0123 in 



c = 0.2H- B/lbm 
P 



First Test 



it„ 



k-1 



surr s k 



) 



535 



.286 



= 389R 



iv 2 = 223.8 V.2H (535*-389) = 1325 ft/sec, 
F_ * 8.5 sin 5 = 0.7U1 lbf 



As the pressure ratio Is = .328 <.53 the mass flow is (see foot note) 

m = . 532 ( r .7)(.0123) s 0 . 01 26g| 
535 sec 



Then for the sub critical pressure ratio, with p* = p. (^3) = 23.7 psia the velocity 



found from eq. (1U-9) is 



v * 

v 2 



• F-ACp^p, ) j Q,.7m-0.0123 (tm.7-23.7) / 



> .0l26 v 
K 32. 2 } 



Second Test 



v 2 = 123H ft/sec , 



n 



i'T 



2 
iv. 



y 0 /iv« = 123U/1325 = .93 



• 535 <M> 



.286 



= 367R 

223.8 .2U (535-367) = 1^19 ft/sec. 
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F_ = 8.5 (sin 6) ,? 0.946 lbf 
R 
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FIG. 14-3 IMPULSE EXPERIMENT 
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Again, the pleasure ratio if th4n^drlH^o*Xt-^ 1 It\ft 



The maaa flow ia (see foot note) 



. J 9 \ m ^: 7} ^ 0123) . • 0,0155 

t sec 



;*35 

Then. with p* ■ 54.7(.53)' » 29.0 psia 



b 0.9U6 - 0,0123^.7-29,0) ; 1309 < ft/9eo> 

2 , . 0155 y . 

1309 n „„ • 

k n " TOT 3 °- 92 ' 



fr'^^i ■ r 

With additional taat data these results can be plotted as a function 6f tha jiozzle 

pressure, Pj»"* or id**^ axit valocity i^v^. "^^^^ x 

Ona poasibla modification of this, akparimarit can ba dona by constructing no*«lea of * 

different lengths or shapes. Than velocity coefficients can ba compared for the different 
nozzle shapes 

A simpler form^of this experiment can be done using .water as the fluid. 
Experiment 1UA 

Measur ement of impure in a Jet of Fluid 1 

1 ■ - 1 «r" 

The object of this experiment is to measure the impul&fe force 
When a jet strikes surfaces of different geometry. 

An extension of the previous experiment can be made by measuring the impulse in the 
jet from the nozzle of Experiment 1H. A system was built, as shown in Fig. 1^3,. The 
nozzle was fastened "to a support, with the jet pointing vertically downward. Placed 
several inches below the nozzle was a balance at which the jet was directed. This per- 
mitted the impulsive forfce (reaction force in the jet) to be measured by means of the 
balance. The objective of the experiment was to measure the force with several shapes used 
to change the impulse. Baa^d on measurements of "T^ p 1 and p ? the mass and velocity 

(otf momentum) in the jet could be determined. These Are idealised values. Then the measured 

-<- 

impulse force could be compared to the value computed from the pressure ^nd temperature 
measurements. Fiv$ shapes were used. These are shown in Fig. 1U-4. They vary from a simple 
flat plate to surfaces with two dimensional flow dividers. 

f. 

Method of calculation 

For each test, the pressure in the nozzle, p t was set, and readings 6f P 8Urr 

and the force on the balance were recorded. 

To compute the theoretical force, the momentum of the jet at the nozzle exit was found 
as the product f of mass, A, times velocity. , 
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H«r«, r«f«rring to Exp«rin»Tit It, «q. (!«♦-«»), (1U-5), (l»»-7) and the g«» law. 
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For casrffc in whifch the pressure ratio is lean than critical determine Ae mass 

flow as ~~ /* ! -' - ■■ ■■ - - ----- - 

p A 

ifi = 0.532 



The forc< depends on the geometry of the object placed on the balance. For steady 
flow, the force can be found by applying the momentum equation to a cont ^ol volAne. As 
the jet is not confined the pressure forces are balanced, thus the force against the object 



is 



F = momentum flux in vertical direciton flowing out 
- momentum flux in vertical direction flowing in 



For the flat plate, case A, the efflux of momentum is; horizontal so ifcv^ = 0 , . 
and % 

F = mv 2 - M-10 

where is the nozzle Velocity or gas velocity in y- 

direction entering the control volume. 4 
This is also true for configuration B. The sphere andsplit ^hrough, ' cgses C and D' } 
have y- momentum flowing ou£ of the control volume. It. is 

momentum out = ~^ v 2 

Consequently, as" the in-flux of momentum is ^ v 2* tota * ^ orce * s » 

F = 2mv 2 , W-ll 

With the exit reduced to 50°, configuration £, the efflux was 

mv 2 sin 50° = mv 2 (.7660) ' 



and the force on the balance- 



F = mv 2 + mv 2 ( ,7660) 
. = 1.7660 mv A 



W-12 



With the several impulse surfaces mounted, ori the balance, the actqal impulse force 

F* " can be measured as a function of the pressure in the nozzle. From the equations for 

R . 

ideal^flow, the ideal nozzle exit velocity, iy^, and mass flow can be found as a function 

of nozzle .pressure, and the ideal impulse computed, F. From these commutations a ratio of 

measured force F D to ideal force F, can be used to express the effectiveness of each 

R ^ 
surface. Then thede ratios can^be plotted as a function of the ideal exit velocity, iv . 
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Sampl? results 

experiments were run it&JLng t;lff<f 0 . 18$ in ^(Hamater nozzle , and the fivti ehapos" ahqyn „in 
Fig. 1M-H. As the tpmp©r«lur4/^>{.JtKe rtd^il«.i^«Ine a ^QoaiTly constant it was possible to ; 
generate curvea pf i^eal ^xit : ;y<i^ as^a function of nozzle pressure and 

surface shape. Curves of itipu.ttfflt forcjfc f dta > the - if i*?e shapes and ideal exit y^locjLty for the: 

nozzle are shown in Fig. JL^S* fcfc each- experiment , the ;ldeal impulse- force was taken for 

: : ?"v • ■ i? V - ' ' ^ 

the measured nozzle pressure from Jthe appropriate curve of fig. m-5, This was then compared 

to the measured force., .in the r&t£?>. Fp/F. .tfsing'the *id£al velocity taken from the curve in 

Fig. 1M-5, this ratio was plotted .to produce* the points in Fig. . % / ' 

These data are from a student "conducted- experiment t vlnd no attempt is made here to , 

evaluate the results. However*, . there is- cohsfistencjy ftmong data for each surface, and it is 

evident in every set that there , Is an I^cf^ease' tih v '%tf*&iwxiWB : ^av lv 2 lUOtf ;Pt/sfe<r^ 

Consequently* this inflection was incorporated v in the M ^|pd" curves drawn thr6ugh. the data. 

It is, more pronounced for cases A an&*> B, 

Experiment mB ^ s - u # - ' <■ ~. 

Distribution of a Fyee, Jet in Air y mm ; - - - 

The object is to measure the shape of the plume of l% 
a jei of air into air, * 

Using the nozzle system of Experiment l l * , it i& . possible to map the distribution ( of 
the jet in a inborn of still air. Vfith a simple hot-wire system , velocity can be measured in 
planes aligned normal t6 the nozzle axis and at different distafices -from the nbzzle exit.. 
This provides a thre<? dimensibnal distribution o£ velocity in the plume p£ the J^t . ; ljaps f 
showing lines of constant velocity can be •Constructed, ' t ',. ' 

' ■ A .'simpler version of the above \ chxV be done by . constructing a square about m 4 ft ^ .Hft> M 
of lig^ht piping, mounted 'oft - a portable stand. Twine is. thej^attachred 'across the piping to. 



form a square grid as shown in Fig. m~7. At the intersection of the cro$s .strings , small, 
tufts of light y$rn are fastened; When the air jet is started .some -of the tuf t& ,on. the frame 
(l will be affected by the air movement;;.- By moving the frame along/the Itxis of the 1 jet', cf< 

qualitative distribution of th$ jet plume can be constructed.,^ ♦ ' ' 

* . - • * . '* ' " . -.'>' * ' - 

In each of the cases, the nozzle pressure can be changed aitd the jet di^tixibut iort .J. ; 

mapped as a functipn 6f the nozzle pressure. . # ^ ' ^ * ' 

, _ ■ ' , : r*> ;> h ■■ • . * ^ . f .. 

Experiment ^1MC , " ■ 4 . • "'^ 

Noise , in a Jet ' - ' 1 * ' ' ' ' ' ' - " / 

* • " ■ * • 1 ■ ' 

1 The objective is to measure thg^soun^ level -and * • . •; * <* 

< ^ distributiop of sound level In *a jet plujfie,. . / ' ^ * 



Another experiment was run with the nozzl^ to determine the .dW trit»uti(5n pf souftd 
The nozzle Was run in free airland tAVe boom probed witlx a ainall 1 micj^tfpb'one and a .^'odnd* # * 
level meter.* From the readings maps were prepared to'sUpw tlie leyel'^f sountl lh..J;he j<^t, m 

and in' the roQra. Th© distribution of sound Vias then compared^ to 'thfe distribution oF Jiifc ~" 

4 . " \ - • .* ^> ••" . 

vel6city. A line of maximum .sdond intensity wa?, found starting at.-.the nozzle- c^riter.lift^ VJ. 



•tid spreading conical,ly, until at a diatanoe 2.«* ft. fr^ the exit plan, the pone radius ( •? i 

was 9 feat. Sound level map's ware prepared in two dimensions for a no»zle pressure of « A* 

- - - m 



• In addition, two mufflers were constructed, applying 2 inches of fibreglass insulation - 
to wire cloth. One muff ler was cylindrical, 3 feet long and 22 inches inside diameter. The 

eeoond wa* donie , 33 inches long^with an inside -diameter of H-l/2 inches on one end and - 

19-1/2 inches on the other. Sound suppression was disappointing, and the velocity patterns 
were 'altered by both suppressors. At the exit of the conic muffler sound levels changed, 
periodically * It was also observed that the flow as observed by a hot-wit* trace on an 
osciliiscope changed from laminar to turbulent to laminar as" the sound level Ranged/ * 



Experiment 1UD _ 

A Simple Shadowgraph 



The objective is to view a compressible jet 
using a simple, shadowgraph. 



■ When a gas flows at a high speed, the density variations in the gas become sufficiently 
large to affect the index of refraction. In fact density variation due to Wiperature dif- 
ferences cause changes in the index of refraction so large that nearly everyone has observed 
them. These are the distortion waves seen rising from V road on a hot day, tor the waves 
that can be. Observed rising from the shadow of your hand in sunlight. This is the principle 
of the shadowgraph, observing the distortion as shadows cast by a .source of parallel light. 
There are more sophisticated optical systems, for example the schlieren system or the 

interferometer,. .. ^ 

A simple shadowgraph can be constructed to observe the wave patterns In a supersonifc 
jet using a slide projector as a light source. v Place the projector about tifh to twenty 
$ feet from the nozzle, jet with the light shining through the jet at right angles to it. Then 
place a ground glass or white piper, close to the jet and in the light path. It may be neces- 
sary to put the paper 1/2 inch from the jet- center line to get the shadow iu^ focus . Now run 
the nozzle with nozzle pressures more than 1.9 times atmospheric to insure supersonic flow. 
A diamond pattern of shock and. expansion waves will be observed in the shadow cast on the 
ground glass. Because the nozzle is small, the shadowgraph is small and faint, but it can 
be observed by magnifying it; photographed with a close-up lfens ; or displayed on a television 
monitor with close-up lens in the, video camera. ^ p 
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